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ABSTRACT 


The performance of mechanically fastened composite joints was 
studied. Specifically, a single-bolt connector was modeled as a pin- 
loaded, infinite plate. The model that was developed used two 
dimensional complex variable elasticity techniques combined with a 
boundary collocation procedure to produce solutions for the problem. 
Through iteration the boundary conditions were satisfied and the 
stresses in the plate were calculated. Several graphite epoxy laminates 
were studied. In addition, parameters such as the pin modulus, 
coefficient of friction, and pin/plate clearance were varied. 

Conclusions drawn from this study indicate: 1) The material 

properties (i.e., laminate configuration) of the plate alter the stress 
state and, for highly orthotropic materials, the contact stress deviates 
greatly from the cosinusoidal distribution often assumed; 2) Friction 
plays a major role in the distribution of stresses in the plate; 3) 
Reversing the load direction also greatly effects the stress 
distribution in the plate; 4) Clearance (or interference) fits change 
the contact angle and thus the location of the peak hoop stress; 5) A 
rigid pin appears to be a good assumption for typical material systems. 
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Chapter 1 
INTRODUCTION 

Background 

The joining of two composite laminates In structural applications 
can be accomplished using several methods. Adhesive bonding and bolted 
joints are by far the most common. While adhesive bonds are simple to 
fabricate, they have been found to be weak in shear [1] and have the 
further disadvantage that they cannot be readily disassembled. Mechani- 
cal (bolted) joints, on the other hand, are easy to disassemble but can 
lead to undesirable stress concentrations. Since mechanical joints are 
the only form of joining that permit disassembly , they have received 
much consideration in previous research work. 

Figure la shows a simple, double-lap joint bolted joint. In this 
type of joint the load is transferred from the central plate, or inner 
lap, through a bolt to two support plates, or outer laps. Herein the 
total load transfered by the joint is denoted by P and it is assumed the 
two support plates each react one-half the load, or F*/2 . 

The general areas of interest for bolted joints are: determining 

the stress distribution in the plates; determining the failure load; and 
determining the faHure mode. Of course the answers to these questions 
depend on the geometry and material properties of the joint. Experimen- 
tal techniques are quite costly and so various analytical techniques 
have been proposed in an attempt to answer these questions. To simplify 
the analysis, it is often assumed that a pin-loaded plate accurately 
represents the central plate of the joint. By symmetry of the joint, it 
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Fig. 1 Mechanical Joint 

a ) Double-Lap Joint 

b ) Pin-Loaded Plate 
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is generally assumed an analysis of the central plate provides answers 
for the support plates. Figure lb shows the central plate isolated from 
the joint. The bolt has been replaced by a simple pin. It is assumed 
the plate is loaded by a known pin displacement and the effects of this 
displacement are reacted by a far-field load P. This is the model used 
by several previous investigators and in the work presented here. 
However, even though the overall model used here is similar to that of 
other investigators, the details of the model to be discussed are unique 
and felt to represent a step forward in analysis. To substantiate this 
claim, and to provide a general framework for the study, it is prudent 
to review previous resparch involving bolted joints. 

The analysis of bolted joints has been conducted both analytically 
and experimentally by many researchers. Analytical investigations have 
included plane-stress elasticity solutions as well as two- and three- 
dimensional finite-element analyses. The elasticity solutions generally 
assume a pinned connection rather than a bolted connection. This is due 
to the two dimensional limitation of the elasticity solutions. Of 
course with this assumption, any three-dimensional or through-the-thick- 
ness effects are not accounted for. Therefore, bolt clamping force and 
interlaminar effects in composites, for example, are not accounted for 
with any of these elasticity solutions. Recently three dimensional 
finite elements have been used in an attempt to account for the through- 
the-thickness effects of bolted joints. Tne following is a summary of 
some of the work done in the area of bolted joints. 


Literature Review 


One of the first analytical studies was done by Bickley in 1928 
[2], Using an isotropic elasticity approach, Bickley approximated the 
contact pressure caused by the pin with a sinusoidal traction applied 
over one-half the boundary of the hole. This approximation has been 
proven -ecentiy to be a good approximation for isotropic plates if fric- 
tion In the contact region is neglected as shown by de Jong [3j. 

In an early effort to intioduce finite-width effects, Theocaris [4] 
used an elasticity solution to model a rigid pin in contact with an 
isotropic plate. More recently Oplinger and Gandhi [5] used a least- 
squares boundary collocation technique to introduce finite geometry and 
multiple-hole configurations into their analysis of pin loaded orthotro- 
pic plates. Superposition of two infinite plate solutions was used by 
de Jong [3] to approximate the finite geometry effects of the orthotro- 
pic plates in his study. In a later report [6], de Jong used an infi- 
nite plate solution to study the effect that load direction has on 
orthotropic plates. When the loading direction dia not coincide with 
one of the principle material directions, the shear-coupling effects 
were found to alter the stress distributions in the plates. 

All of the analyses discussed so far neglected friction in tne 
contact region. The friction between the pin and plate is a non-con- 
servative force and is therefore difficult to include in any analysis. 

In one of their studies Oplinger and Gandhi [7] introduced friction into 
the contact regions. They produced preliminary results of the limiting 
cases of slip in the entire contact region and no-slip in the entire 
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region. Oplinger and Gandhi did find that the presence of friction 
altered the predicted failure locations in some of the laminates they 
studi ed. 

Elsewhere, researchers at the Indian Institute of Science have 
introduced variables such as pin elasticity, friction, and Interference 
(or clearance) fits into the problem for isotropic plates. In the 
findings of two reports [8, 9], the researchers were ultimately able to 
analyze an elastic pin inserted in an isotropic plate. The analysis 
included an interference fit and finite friction in the contact 
region. The analysis was taken only as far as the onset of separation 
for the interference fit. A large amount of information was compiled by 
the researchers and was presented in the form of nondimensional ized 
plots. In a third report [10], the researchers at the Indian Institute 
of Science studied the stresses in an orthotropic plate loaded by a 
rigid, frictionless pin. In this work the researchers focused on the 
effects of orthotropic elasticity rather than the details of pin/hole 
interaction. 

The most recent elasticity approach including friction in tne 
contact region was published by de Jong [11]. In a modification of his 
earlier work, de Jong Included friction in the slip zone but was forced 
to use an empirical curve-smoothing procedure In a region which he 
defined as the "release area". In defining the release area, de Jong 
relaxed the tangential displacement boundary condition in the no-slip 
region and applied a smooth, anti symmetric shear stress which was a 
function of material properties and loading direction. With this analy- 
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sis de Jong showed that friction had a significant Influence on the 
stresses In a pin-loadea composite panel. 

The papers discussed so far that have included friction have done 
so only for monotonic loading. Reversal of the load direction (e.g., 
fatigue loading) reveals the nonconservative nature of friction. This 
complicates the problem significantly. In fact, many researchers have 
chosen to ignore friction entirely in an effort to make the problem 
mathematically tractible. In particular, finite-element models are 
often used without including friction. 

In an early orthotropic, two-dimensional finite-element solution, 
Waszczak and Cruse [12] used an assumed traction boundary condition 
similar to that chosen by Bickley [2], Six laminates were analyzed and 
a distortional energy failure criterion was used to predict failure. 

The researchers stated that they obtained good agreement with experimen- 
tal results for the prediction of the failure mode but they found gen- 
erally poor agreement for the failure load. This study accounted for 
finite-width and end effects but did not accurately model the tractions 
In the contact region. 

Agarwal [13] used a displacement boundary condition to model a 
rigid, frictionless, pin-loaded plate. The results of the finite ele- 
ment model were used in conjunction with an average stress failure 
criterion to predict the failure of single fastener joints. Agreement 
between analytical and experimental results was found to be better than 
Waszczak and Cruse [12]. Exceptions to this were the [+45] s and [0/90] $ 
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laminates. With these laminates the predicted strengths were approxi- 
mately 50% of the experimentally measured values. 

The finite element model by Chang et al . [14] used a traction 
boundary condition similar to that of Waszczak and Cruse [12], A modi- 
fied average stress failure criterion was used to predict the failure of 
laminates tested by other researchers [12, 13, 15]. The authors noted 
that all of the predictions were conservative for the [±45] $ and [0/90] s 
laminates. In papers [12-14] the stress distributions were not compared 
with previous solutions. In fact, two of the papers [12, 14] have used 
traction boundary conditions which other researchers have disputed as 
being valid for orthotropic materials. This may help explain the 
incorrect failure predictions of ref. 12 and 14 for the highly orthotro- 
pic [+45] $ and [0/90] s laminates. Despite this problem, improvements on 
the failure hypothesis seem to have been made. 

Crews [16] modeled a flexible pin in an orthotropic plate using a 
finite element technique. By using a flexible pin the frictionless 
displacement boundary conditions were satisfied. In general, the 
results that Crews obtained agreed with those of others [5, 8], Crews' 
paper primarily studied the effects of finite geometry on the stress 
distributions in single pin connectors. No attempt was made to predict 
failure. 

Including friction in any finite element model is difficult and 
often requires a reformulation of the energy equations for the system to 
be analyzed. It may be possible to avoid this by using existing finite 
element techniques which satisfy the boundary conditions indirectly or 
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approximately. There are few finite element solutions which Include 
friction. However In his Ph.D. dissertation, Wilkinson [17] Included 
friction In his finite element model. He did not compare his solutions 
with any previously published solutions and many of nls results seem to 
be Inconsistent with the findings of other researchers. For example, in 
a paper published by Rowlands, Rahman, Wilkinson, and Chiang [18], the 
authors stated that the radial stress component in the contact region is 
''nly slightly changed with large changes in the coefficient of fric- 
tion. This contradicts the findings of earlier researchers [7, 11], 

Recently, Matthews, Wor.g and Chryssafitis [19] developed an isopar- 
ametric brick element to analyze the effect of bolt clamping pressure 
and other through-the-thlckness effects. Using only one element through 
the thickness the authors were unable to account for interlaminar 
effects. A rigid pin with no friction was used to apply a load to the 
orthotropic plate. Matthews et al . found that bolt clamping pressure 
altered the stress distribution in the plate. 

Several other researchers [20-23] have used finite element methods 
to study bolted joints. However, no one, using any method, has success- 
fully treated the case of an elastic pin in an orthotropic plate with 
finite friction in the contact region. In addition, the effects of 
pin/hole tolerance have not been addressed for the case of an orthotro- 
pic plate. Most researchers have assumed a push-fit pin. 

Where analytical methods have failed, experimental investigations 
have been conducted in an 2 ttempt to obtain information. Among the 
first attempts to experimentally determine the stresses in pin-loaded, 
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isotropic plates were the studies undertaken independently by Frocht and 
Hill. Their results were combined In a single paper appearing In 1940 
[24] and later in Frocht's book on prtotoelastlclty [25], Frocht chose 
to analyze the problem using photoelastic methods while Hill used strain 
gages. The strain gages used had relatively large gage lengths (0.5" 
and 1.0") and Interpolation was required to find the strains at the edge 
of the hole. Several diameter-to-width ratios and clearances were 
tested by both authors and, in general, good agreement between both 
investigations was found. 

Nisida and Saito [26] developed an interferometric method which 
they used to predict the stresses in a pin-loaded isotropic plate. They 
presented limited results on the stresses in the plate seemingly to 
prove that tneir interferometric method was valid and could be used as a 
future research tool. They did rote, however, that the maximum hoop 
stress occurred at an angle about +85° from the loading direction. This 
seems to confirm the results of many analytical investigations. 

Wilkinson, Fuchs, and Rowlands [27] used photomechanical (moire ard 
holography) techniques to obtain the response of pin-loaded holes in 
Sitka spruce and Balsa wood. The results were compared with previous 
analytical work [16, 17] at a li-ited number of locations in the 
plate. For the locations that were compared, the agreement between the 
analytical and experimental studies was found to be good. 

Recently Prabhakaran [28] and Hyer and Liu [29] have used orthotro- 
pic photoelastic techniques to study the stress distribution In glass- 
fiber reinforced plastics. While the subject of orthotropic photoelas- 
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1 1 cl ty Is relatively new, Prabhakaran and H.yer and Liu are using the 
technique to obtain results for the case of pin-loaded holes. The 
results should prove Invaluable in the verification of the analytical 
studies discussed earlier. 

Several other experimental investigations have studied other 
aspects of mechanically fastened joints. Lambert and Brailey investiga- 
ted the importance of friction between the pin and plate [30], In their 
experimental analysis, Lambert and Brailey Investigated the effects of 
slip and separation in the interference fit joints of isotropic 
plates. They found that a coefficient of friction equal to 0.3 existed 
for the Araldite specimens they used. Alterations in the coefficient of 
friction were obtained by coating the surfaces with a lubricant. It was 
found that changing the coefficient of friction significantly changed 
the shear stress concentration factor. 

In studies involving through-the-thickness effects, Quinn and 
Matthews varied the stacking sequence of 8 ply, quasi -isotropic, glass- 
epoxy specimens [31]. Eight stacking sequences were tested and the 
results showed that ultimate failure was dependent on the stacking 
sequence. The results seemed to suggest that placing the 90° plys near 
the surface of the panel Increased the bearing strength of the pinned 
joints. 

Stockdale and Matthews studied the effect of clamping pressure 
resulting from torque applied to the bolt in a joint [32]. Two 0/90 
laminate configurations were tested with various washer sizes an' bolt 
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torques. Stockdale and Matthews found that large washers and high bolt 
torques led to improvements In the strength of the joints. 

In a related study, Shivakumar and Crews developed an equation to 
describe the relaxation that occurs in resin-matrix composites where 
bolt torque is a factor [33]. The authors exposed several 32 ply, 
graphite-epoxy, bolted joints to various environments for a period of 
one year. Relaxation of the joints was characteri zed using a visco- 
elastic approach. The authors found moderate correlation between exper- 
iment and their viscoelastic theory. 

To obtain higher performance joints, Eisenmann and Leonhardt [34] 
have investigated laminate tailoring concepts. By adding more 45 degree 
plys in the bearing region the authors were able to increase the perfor- 
mance of full scale composite components by as much as 62%. 

Other investigators [35-40] have studied various multiple bolt 
geometries, performance at elevated temperatures, and other factors that 
effect the performance of mechanical joints in composite materials. 

These tests have generally been conducted to develop a data base with 
the hope of developing failure prediction techniques for design pur- 
poses. 

One of the earliest failure prediction techniques was developed by 
Whitney and Nuismer for composite specimens with holes or notches 
[41]. Their hypothesis was based on the idea of a plastic zone found in 
metals. Whitney and Nuismer assumed that a characteristic dimension 
existed over which the laminate must be critically stressed in order for 
failure to occur. This theory is the average stress failure theory 
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already discussed In the review of the papers by Agarwal [13] and Chang 
et al. [14], In these cases the stress distribution was calculated 
using finite-element models but Garbo and Ogonowski [15] used an elas- 
ticity solution for stress predictions. In an extensive study, Garbo 
and Ogonowski developed an analysis technique which predicted the fail- 
ure of bolted joints In composite materials. The authors used a complex 
variable elasticity approach to predict the stresses and a choice of 
several failure criterion to obtain failure predictions. The elasticity 
solution, however, relied on an assumed boundary traction similar to 
that of Blckley [2]. This boundary traction has been shown to be 
Inaccurate for some orthotropic materials. Both de Jong [5] and Crews 
[16] have shown this. This fact casts some doubt on the technique's 
ability to predict failure accurately. 

Hart-Smith has also developed a capability to predict failure In 
mechanical joints [42]. Hart-Smith conducted many experimental tests 
with single and multi pi e-b._.lt joints. From these tests he developed a 
failure hypothesis based on stress concentration factors. The system 
that Hart-Smith developed relied heavily on a large data base and aocs 
not appear to have the flexibility required to handle joint configura- 
tions other than those already tested. 

Using a somewhat different approach, Col lings [43] has developed a 
failure prediction ~apab111ty based on his experience with bolted 
joints. Col 1 lugs chose a semi -empl rl cal approach which used Bickley's 
loading on a hole [2] to calculate the contact pressure. However, 
rather than Ignoring friction, Collings used a sinusoidal friction 
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distribution in the contact area. Thus, Col lings' calculation of 
stresses would seem to be a more accurate representation of actual 
bolted joints than the representation of Garbo and Ogonowskl who Ignored 
friction. The reasoning for some of Col lings empirical deductions seem 
to be based on a strength of materials approach that relies heavily on 
his experience In the area of bolted joints. Using tensile and com- 
pressive strength data obtained from coupon tests, Col lings shows good 
agreement between his failure predictions and experimental results. 

All of the failure predictions just discussed required some degree 
of empirical formulation and they lack the flexibility which Is required 
for general design purposes. The methods do not use accurate descrip- 
tions of the stress state in bolted composite laminates. Rather, they 
make assumptions which are enforced independent of the particular joint 
geometry or material. Since an accurate description of stresses Is the 
starting point of all failure analyses, It seems that a method dependent 
on the specific parameters of the joint would be more appropriate. In 
addition to the accurate prediction of stresses, it seems that a more 
flexible failure hypothesis needs to be developed. This has been 
pointed out in other reviews of past research. Recent reviews by Godwin 
and Matthews [44] as well as Opllnger [45] have been written in the area 
of bolted composite joints. Godwin and Matthews suggest that a factor 
of safety of two must be used In current designs and the review by 
Opllnger demonstrates the complexity of current failure prediction 
techniques. These reviews indicate that more work is required before a 


comprehensive understanding of the performance of composite bolted 
joints Is to be achieved. 

Objective 

The failure prediction techniques developed to date have lacked the 
flexibility required for design purposes and, in many cases, they do not 
predict the failure accurately (e.g., the [0/90] $ and [+45] s laminates 
of refs. 12 and 14). For accurate failure predictions for bolted 
joints, two steps are necessary. First it is necessary to accurately 
predict the stresses in the joint. Second, it is necessary to use a 
valid failure hypothesis. Without the first step it is unl^xely that 
accurate failure predictions can be made. This thesis reports the 
development of a more accurate method for predicting the stress distri- 
bution near holes in composite joints. 

From the literature review it is evident that various boundary 
conditions at the hole edge have been used to predict the stresses in 
composite joints. This is true for both the finite element approaches 
and the elasticity techniques. Often these boundary conditions are 
unjustified and may in fact be in error. Parameters such as clearance, 
pin flexibility, and friction are often ignored. Furthermore, there has 
been no analysis which includes all of these parameters for composite 
joints. A two-dimensional analysis of pin-loaded joints can be used to 
study the effects of these parameters. This avoids introducing the 
complexity of a three-dimensional analysis into the problem, yet new 
results regarding these parameters can be obtained. This is the 


approach taken In the study here. In particular, It will be assumed the 
pin Is in a state of plane strain and the plate Is In a state of gen- 
eralized plane stress. The plate Is assumed to be infinite In extent, 
using these assumptions, an analysis, based on complex variable elastic- 
ity, Is developed which can accommodate pin/hole clearance, pin flexi- 
bility, and friction between the pin and the hole. By varying these 
effects Independently, a better understanding of the problem will be 
possible and the Importance of these often-ignored effects can be 
assessed. 

In the succeeding five chapters the problem Is formulated, the 
equations governing the pin, the plate, and their interaction are devel- 
oped, and numerical results are presented. Chapter 2 presents the 
nomenclature and geometry used, and basically defines the problem. The 
method of solution Is outlined. The problem, being a contact elasticity 
problem, belongs to a class of problems for which there are few closed 
form solutions. Numerical techniques are usually required If quantita- 
tive Information Is desired. Here no attempts were made to find closed- 
form solutions, rather It was assumed from the start that numerical 
techniques would be necessary. Chapter 2 outlines the numerical scheme 
used here, a scheme based on the collocation method and an Iteration 
process. The elasticity solution for the pin is derived In Chapter 3 
and Chapter 4 presents the details of the elasticity solution for the 
orthotropic plate. That a separate analysis Is needed for each compo- 
nent will be made dear In Chapter 2. Chapter 5 discusses details of 
the numerical method, Including convergence of botn the collocation 
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method and the Iteration procedure. Finally, Chapter 6 presents some 
very Interesting results regarding the stress distributions around the 
holes. The effects of pin flexibility, pin/hole clearance, friction, 
and degree of plate orthotropy on the stresses are demonstrated. Chap- 
ter 7 concludes the work and presents recommendations for further 
studies. 


Chapter 2 

PROBLEM DEFINITION AND METHOD OF SOLUTION 
The problem of a pin-loaded connector is a geometrically nonlinear, 
fourth-order boundary value problem with non-conservative forces acting 
between the pin and the plate. Figure 2 illustrates the nonlinear, 
nonconservative nature of the problem. In this figure a [02/+45] $ 
graphite-epoxy laminate is loaded using two different assumptions. In 
the first case the clearance between the pin and plate is zero and 
friction is neglected. The condition of zero clearance is referred to 
as a push fit and is an often-used assumption. This results in a linear 
relationship between the pin load and pin displacement. Due to the lack 
of friction, the load-displacement relationship is independent of load- 
ing direction. This case is compared to the nonlinear and path-depen- 
dent case of a clearance-fit pin and friction. The nonlinearity in the 
load-displacement relation is due to the clearance. With clearance, the 
contact region between the pin and the plate changes as the load 
increases. This changing contact region leads to the geometric non- 
linearity. Friction accounts for the fact that the loading and unload- 
ing paths are different. The direction of the frictional forces depends 
on whether the connector is being loaded or whether it is being 
unloaded. Fortunately, for any given value of pin displacement, the 
problem is linear elastic. Through the use of the principle of super- 
position, a collocation procedure, and iteration, a solution to tne 
problem can be found. 


17 


Pin 

Load 

(lbs.) 



Fig. 2 Illustration of 
Nature of the P 




19 


Plate and pin models 

In modeling the plate, a two-dimensional analysis Is easily justi- 
fied through the assumptions of generalized plane stress theory. Gen- 
eralized plane stress Is consistent with classical lamination theory and 
therefore can be readily applied to composite materials. Of course this 
will mask any Interlaminar effects In composite materials. However, as 
stated earlier, the effects of pin elasticity, clearance, and friction 
can be studied using two-dimensional models and the inclusion of these 
effects is a step forward. The use of an Infinite plate greatly simpli- 
fies the elasticity analysis. Furthermore, meaningful parametric 
studies can be made regarding the effects of pin elasticity, friction, 
and clearance even for the infinite plate. It has been shown that for 
plate-width to hole-diameter ratios exceeding 8 or 10 and end-distance 
to hole-diameter ratios exceeding 4, the effects of finite geometry are 
negl igible. 

An accurate two-dimensional model of the pin is more difficult to 
justify. A major concern is that the pin must be in equilibrium. The 
reactive forces on the pin, from the plate, represent a nonequilibrating 
system on the pin boundary. Figure 3a shows a free-body diagram for a 
typical pin in a double lap joint. In this figure the load from the 
plate, P, is equilibrated by the forces from the two support laps. To 
reduce the pin to a two-dimensional problem, a disc the thickness of the 
Inner lap is removed from the pin. A free body diagram showing only the 
forces acting in the plane of the disc is given in fig. 3b. It is 
apparent that the load from the plate is equilibrated by the shear 
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Fig, 3 Pin Equilibrium 
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stresses acting on the upper and lower surfaces of the disc. The dis- 
tribution of these shear stresses is not known aid only a three-dimen- 
sional analysis can yield this information. As a first approximation, 
it is assumed that the shear is uniformly distributed over the upper and 
lower surfaces of the disc. As the thickness of the disc Is decreased, 
the stress distribution due to these shear stresses approaches the 
stress distribution caused by a uniform body force within the disc. The 
pin model adopted in this study therefore assumes that the tractions on 
the boundary of the pin are equilibrated by a constant body force within 
the pin. Thus, it seems that the accuracy of this model is limited only 
by the assumption that the shear stresses vary uniformly over the upper 
and lower surfaces of the disc. By using this approximation, however, 
the problem is reduced significantly and only the boundary tractions 
which generate the load P are left as unknowns. 

Coordinate system and nomenclature 

The coordinate systems used in this analysis are shown In fig. 4. 
The x-y axes coincide with the principal material axes of the orthotro- 
pic plate. In this analysis, the load Is introduced into the plate by a 
rigid body displacement of the pin in the positive x direction. 

The r-0 axes are typical cylindrical coordinates, positive 9 being 
measured counterclockwise from the +x axis. The boundary conditions at 
the hole edge and at the pin boundary will be written In the polar 
coordinate system. Numerical results will be presented using the polar 
coordinate system. The displacement in the x direction is denoted by u 
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while the displacement In the y direction Is denoted by v. The radial 
displacement, which is assumed positive In the +r direction, Is denoted 
by u r . The tangential displacement, which Is positive In the +9 direc- 
tion, Is denoted by u 0 . The three components of stress In the x-y 
and r-0 systems have the usual definitions and are denoted as 

o , o , x , and o , a nn , and x respectively. The strains In the 
xx yy xy rr 09 re 

two systems are: e xx> e yy , y xy and e rp , e 00 , y rQ . If the circle In 
fig. 4 represents the pin, then the boundary tractions act on the out- 
side of the circle. If the circle represents the hole In the plate, 
then the boundary tractions act on the Inside of the circle. Although 
this may be obvious and is the convention dictated by elasticity, there 
can be somp confusion on this point when considering traction and dis- 
placement conditions at the contact region between the pin and the 
plate. 

Figure 5 defines the pin/hole clearance, \, and the rigid body pin 
displacement, 6. The clearance parameter is positive for a clearance 
fit, zero for a push fit, and negative for an Interference fit. As will 
be seen, the push fit case, \ = 0, Is a special case. Both the loading 
and unloading cases require a positive pin displacement 6. The case of 
loading assumes that the pin displacement has been increased from some 
lower value of 6, while the case of unloading assumes that the pin 
displacement has been decreased from some higher value of 6. For exam- 
ple, loading occurs if the pin displacement increases from 
6 = 0.02 to 6 = 0.03. On the other hand, unloading occurs when the pin 
displacement decreases from 6 = 0.04 to 6 - 0.03. In both cases the 
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problem can be analyzeo for 6 = 0.03. In the methodology here, the only 
difference between the two cases will be the sign on the shear stress at 
the pin/plate interface. 

Conditions between pin and plate 

Figure 6 shows the different regions of the pin/plate Interface 
where various boundary conditions and Interface conditions exist. As 
can be seen in the figure, three important regions exist. The regions 
are denoted as regions I, II and III and these are defined In the fig- 
ure. The actual extent of each of the three regions depends on the 
pin's elastic response, the plate's elastic response, the pin clearance, 
the pin displacement, and the coefficient of friction, p. However, 
independent of the actual extent, rich region has its own unique bound- 
ary and/or interface conditions. These are as follows: Let the radial 

and tangential elastic displacements of the pin and plate be defined, 
respectively, as 

u r » * u r » U Q 

pin pin plate plate 

Then in region I, referred to as the no-slip region, at the boundary: 

u = u +6 cos 9 - \ ( 1 ) 

r plate r pin 

- a < 9 < a 

u = u - 6 sin 9 . 

plate pin 


( 2 ) 
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In region II, where there Is slip: 

u = u +6 cos 0 

r plate r pin 



a < e < p 
and 

- 0 < 0 < - a 


( 3 ) 

( 4 ) 


In region III, the no-contact region: 


o = 0 
rr 



9 < 0 < 2n - p 


( 5 ) 

( 6 ) 


Equations 1 and 3 express continuity of radial displacements in the 
contact region and eq. 2 expresses continuity of tangential displacement 
in the no-slip region. Equation 4 is the Coulomb friction law assumed 
to hold In the slip region. Of course eqs. 5 and 6 express the trac- 
tion- free conditions of the no-contact region. The second and third 
tennis on the right side of eq. 1, and the second term on the right side 
of eq. 2 represent the effects of rigid boay pin displacement and the 
effects of clearance. The coefficient of friction appears in region 
II. Therefore, the loading direction determines the sign of the shear 

stress in this region. For loading, x- Is positive (p is negative) and 

rtf 

for unloading, x r0 is negative (p is positive). 

The coefficient of friction, p, Is often assumed to be only posi- 
tive. However, for this report, the coefficient can change signs. This 
is merely a simple way of accounting for the loading/unloading cases. 


Method of solution 


\ 
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Referring to fig. 6, the unknown boundary tractions act over 
regions I and II. By equilibrium considerations, the boundary tractions 
on the pin are identical to the boundary tractions acting on the hole in 
the plate. The problem reduces to finding these tractions. Once the 
tractions on the hole boundary are known, the stresses in the entire 
domain of the plate can be determined. If they are of interest, the 
same can be said of the stresses within the pin. Of course the elas- 
ticity of the pin, the elasticity of the plate, the pin/hole clearance, 
the pin displacement, and the level of friction, all couple to determine 
these tractions, and, as mentioned previously, the actual domains of 
regions I, II, and III. How can these tractions, which are the key to 
the problem, be found? The answer lies in the use of superposition. 

To make use of superposition, the unknown tractions on the boundary 
of the pin and plate are taken to be represented in the form of a com- 
plex Fourier series with unknown coefficients A^. The series is: 

< N - iT > pin ■ < N -' T > plate = T\ e ,ke (7) 

where N = the normal traction = o rr @ the boundary 

T = the tangential traction = T f0 @ the boundary 

1 = v^T 

A k are the unknown coefficients of the Fourier series. 
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There are several points regarding eq. 7 that should be discussed. 

First, the equation Indicates normal and tangential tractions on the pin 
edge are equal to the normal and tangential tractions on the hole 
edge. Thus the same Fourier series Is the driving force behind the 
analysis of the Isotropic pin and the distinctly different analysis of 
the orthotropic plate. This Is one of several couplings between the 
plate response and the pin response. Second, since the tangential 
traction must be an odd function of 0 and the radial traction must be an 
even function of 0, the unkno*<n coefficients A k are real. Finally, the 
tractions around the entire circle are represented by continuous func- 
tions of 0. Therefore the coefficients must be such that H and T are 
zero in the domain of 0 corresponding to region III in f 1g. 6. 

1 k Q 

Since they are a legitimate boundary traction, each term e In 

the Fourier series produces a unique set of stresses and displacements 

in the pin and a unique set of stresses and displacements In the 

plate. Through the elasticity techniques presented In the following 

chapters, the stresses and displacements In the pin and plate due to 
1 k 0 

traction e can be found in closed form. If the boundary traction Is 
1 k 0 

A k e , then the components of displacement and stress will be weighted 
by the factor A k - If. for example, the radial displacement due to the 

1 k 0 

traction e Is denoted as (u r ) k , then the radial displacement due to 
1 k0 

traction A k e Is given by 


\< U A 


( 8 ) 
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Considering all terms In eq. 7, tne radial displacement will be the 

1 y Q 

weighted sum of the effects of the Individual tractions e . That Is, 


u r I V u r>k ' < 9 > 

k = _OD 

If all A k were known, then u p could be calculated. The same can be said 
of the tangential displacement, namely 


% ■„ ' Wk ■ (10 > 

k = -“ 

1 k 6 

where (uj is the tangential displacement due to traction e and u_ 

9 k 

is the tangential displacement due to the weighted sum. 

To avcld confusion, no distinction has yet been made between the 

plate and the pin. However, in eqs. 9 and 1C the radial and tangential 

i k 6 

responses due to traction e , specifically (u ) and (u.) , are dif- 

0 k 

ferent for the plate than they are for the pin. This is obvious since a 
pin responds quite differently than a plate to a given set of trac- 
tions. Thus the notation 


u 

r pin 


° pin 
and 

u 

r plate 


ao 


l 

\< u r>k . 

(11) 

Z .CD 

pin 


CD 

I 

A k (u e ) k 

(12) 

= —a> 

pin 


CD 

1 

Wk 

(13) 

s .ao 

K • K plate 
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% = ^ MOi, 

plate k=-» * plate 


( 14 ) 


distinguishes between pin and plate responses. 

It is important to note A k is common to the radial and tangential 

responses of the pin and plate. It is also important to point out that 

i k 0 

since they are a result of boundary traction e , 

(O k > (u e ) k » (u ) k and (u ). 

pin pin r plate y K plate 


are known! The 


problem, as it is being posed here, will be reduced to finding A k . 

To continue with the formalism of the problem, the stresses in the 
pin and plate are given by 


and 


a 

rr pin 

II 

7T 

II 

1 C — 1 

Q 

A. (a ). 
k rr 'k . 

° pin 

(15) 

°06 . 

pin 

II 

TV 

II 

i 8 

8 

A * pin 

(16) 

Ve . 

pi n 

8 

8 i 

il 

il 

' k pin 

00 

rr plate 

n 

7T 

II 

i 8 

8 

A (a ). 
k rr k plate 

(18) 

° ee plate 

00 

- 1 
k=-® 

^k (°eo'k 
k 99 k plate 

(19) 

rG plate 

II 

TV 

II 

1 8 

3 

. t • 
plate 

(20) 


The quantities (o ), , ..., (t V are the stresses In the 

fr k Pin re k plate lfc 

pin and plate, respectively, due to traction e anu they also are 
known. In particular, on the boundary of both the pin and the plate, 
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,T re> 


= N - IT 


boundary 


■ I 

k=-« 


A,, e 


ike 


( 21 ) 


To determine the unknown coefficients A^, the second and third aspects 
of the problem solution, namely a collocation technique and iteration, 
are used. These aspects are explained next. 

As might be expected, finite, rather than infinite sums are used in 
the series in eqs. 7-21. Instead of satisfying the boundary and inter- 
face conditions at each point on the circular boundary (i.e., at an 
infinite number of points), the boundary and interface conditions are 
satisfied at a finite number of points. This technique of satisfying 
required conditions at only a finite number of points is a numerical 
technique referred to as a collocation procedure. The finite number of 
points are referred to as collocation points. The number of terms in 
the finite series is then directly related to the number of collocation 
points. The iterative portion of the solution technique involves deter- 
mining the expanse of regions I and II, i.e., determining a ard p. 

They, like the stresses and displacements, are unknowns. However, if 
they are assigned values, then, since the boundary and interface condi- 
tions are specified in each region, a solution to the problem can be 
found. Starting with assumed values for a and p, iteration is used to 
obtain the correct values. The iteration procedure will be discussed 
shortly. First, the collocation procedure is further explained. 

Assuming a and p are known, then by using finite sums in eqs. 11- 
21, and substituting those sums into the boundary/interface conditions 
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for the various regions, eqs. 1-6, the method for determining the A k 's 
becomes evident. Substituting the finite sums into eqs. 1-6 leads to: 
For region I; 


l V U A 

k=-N 


plate 


N 

= l \ (u ) k + 6 cos 0 - \ 
k=-N K r pin 


N 


N 


l 

k = -N 


V U A 


= l A (u ), - 6 sin 0 

plate k=-N 0 pin 


( 22 ) 

(23) 


For region II : 


N 

l 

k = -N 




plate 


N 

- I 

k = -N 


W\ 


+ 6 cos e - \ 


pm 


k J N Ak((Tre)i( ' " 0 • 


(24) 

(25) 


For region III : 


i y<vA * ° 

k=-N r 

J N \ ( Ve>k 

k=-N 

Here it is assumed there are 2N+1 terms In the sums. With a and 0 
assumed to be known, everything in eqs. 22-27 is known except the 
A k 's, The assumed values of a and p determine over what range of 9 eqs. 


(26) 

(27) 


22 and 23 are enforced, over what other range of e ea$. 24 and 25 are 
enforced, and over what range nf e eqs 26 and 27 are enforced. Since 
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by equilibrium the tractions of the pin and plate edges are identical, 
neither the plate nor the pin has been called out specifically In eqs. 
25-27, Strictly speaking. 


la ). . = (a ) , . = A. cos(ke) \ (28) 

' rr'k oiri v rr'plate k ' ' 1 v ' 

I at 
/ boundary 

^rA pin = AA plate = A k S1n ( ke ) •) ( 29 ) 

Equations 25-27 could stated alternatively as 

N 

l A. (sin(k0) - p cos(ke)) = 0 . (30) 

k— N 

N 

I A cos (k0 ) = 0 (31) 

k=-N * 

N 

l A sin(k0) = 0 . (32) 

k=-N * 


Expanding eqs. 22-27 and rearranging slightly leads to 


U u r ) 0 - (“A J A o + {(u r } -l " (u Al ^ A -1 

r 0 plate r 0 pin 0 r 1 plate r 1 pin 1 


{(Ol " ( U A 1 A l + ••• + ((On - (u L } A 

r 1 plate r 1 pin 1 r M plate r pin " 


6 cos e - x , 


(33) 
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plate ' pin' *° + ll “« , -> plate ' ‘“•’-l p,„* ^ * 


{(Oi - (uj, } A + ... > {(u ) - (u ) N } A 

9 1 plate 9 1 pin 1 9 N plate 9 N pin 


- 6 sin 0, 


{(u ) - (u ) } A + {(u ) , - (u ) , } A . 

lv r'o . _ v r'o . 1 o r'-l , _ v r'-l . ' -1 


plate 


pm 


plate 


pin 


{( u r )i “ ( u r )i } A| + ••• + {(u - (u ),, } A 

r 1 plate r 1 pin 1 r " plate r " pin 


6 cos e - \ 


KVe'o ' “( < ’rr ) «> A o + f (T re>-l * ,l( ‘Vr ) -l |A -l + 

(( Ve'l - “<°rr ) l )A l + — * f< T re>N ‘ “' t ’rT ) N> A N = 0 

t(°rr>o ,A o * l< 0 rr)-l> A -l * * — + <<°rr>N! A N * 0 

^ T reV A o * t^re*-l ,A -l + * ••• * ^ T re*Nl A N ° 0 

From eqs. 33-38 it is obvious that there are 2N+1 unknown A^ 


(34) 


(35) 

(36) 

(37) 

(38) 
■s. If 


the equations are satisfied at 2N+1 discrete points around the circular 
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boundary, I.e., at 2N+1 discrete e locations, then a set of 2N+1 linear 
equations can be formed from which to solve for the A^. This is the 
application of the collocation procedure to the problem at hand. Since 
eq. 7 incorporates the proper symmetry of the problem, collocation over 
one-half the hole circumference, 0 < 0 < n, is all that Is necessary. 
Equation 39 illustrates, in operational form, the collocation procedure 
for finding A Q , A^, A +1 , .... A_ N , A +N . In eq. 39 it has been assumed 
that the collocation points are A0 apart around the circumference of the 
circle. 

The numerical problem now is reduced to the finding of a and 0. 

For given values of a and 0 a solution to eq. 39 is found. The solution 
may not be realistic in terms of the stresses and displacements which 
result. For example, it is physically impossible to have tensile stres- 
ses between the pin and the hole. Also, the shear traction on the 
boundary cannot exceed p times the normal traction. In the solution 
procedure here, trial values of a and p are used. The problem is solved 
via eq. 39 and the values of the normal and shear tractions on the 
boundary examined. If any physical conditions are violated, u and/or p 
are adjusted and eq. 39 is resolved with new values of a and/or 0. 

Again the physical reality of the solution is examined. Within a few 
iterations, values of a and p are obtained that satisfy all the condi- 
tions required. 

It is Important to point out that with a push fit pin, \ = 0 and 
each element on the right-hand-side of eq. 39 changes, for example, by a 
factor of two when 6 is doubled. Each A^ simply doubles and the spatial 


Load Functions 



Equation (39) 
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variation of N-1T, with 9, remains unchanged. Thus a and p are 
unaffected by 6. With \ * 0, doubling 6 does not result in a doubling 
of each element of the right-hand-side of eq. 39. Each A k changes by a 
different amount and thus the spatial variation of N-i T varies. This 
translates into 6 directly affecting a and p. The push fit case is 

therefore quite special and is often referred to as the linear case. 

Chapter 5 will discuss details of the numerical procedure. That 

chapter will discuss, among other things, the specific manner in 

which a and p were determined. The next two chapters are devoted to 

finding the response of the isotropic pin and the orthotropic plate to 

i k 0 

the boundary traction e . Chapter 3 examines the pin while Chapter 4 

presents the analysis of the plate. Both solutions use complex variable 

elasticity. Several books exist on this subject [46-49] but, in gen- 
eral, the details of these types of analyses are not as familiar as 
analyses with real variables. Therefore, some of the details of the 
analyses will be discussed. 


Chapter 3 
PIN ANALYSIS 


For the complete solution to an elasticity problem, three sets of 
equations, in addition to the appropriate boundary conditions, must be 
satisfied. The equations are the equilibrium, constitutive, and com- 
patibility equations. The strain-displacement equations can be consid- 
ered as definitions but certainly they must be satisfied. The boundary 
conditions can be known displacement conditions, known stress condi- 
tions, or a combination of both. In the solution of the problem for 
both the pin and the plate, the governing equations can be simplified 
and combined into one governing biharmonic differential equation. The 
dependent variable in this equation is the Airy stress function. The 
boundary conditions for the problem are satisfied through the use of 
various integrated forms of the solutions of the biharmonic equation. 
This simplifies the solution because after the boundary conditions are 
satisfied, the solution automatically satisfies the governing 
differential equation. The key to both the pin and plate analyses here 
is the development of the stress functions in terms of complex 
variables. This is presented below for the pin. 

The pin is assumed to be isotropic and in a state of plane 
strain. The governing equations are therefore reduced to the following: 


equ'. librium: 


do 5 t 
xx xy _ 

dx dy 


dx do 
dx dy 



(40) 
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constitutive : 


xx 


yy 


xy 


1 

r 

i 

r 

i 


[ < 1 - v2 )°xx 

[(l-v 2 )o yy 


v (l +v )°yy] 

V(ltv) 0 x x] 


compatibility: 


2 »2 

b c d e d y 

xx yy. _xy. 

T . 2 


ay 


ox 


dxdy 


strain-displacement 


(41) 


(42) 


XX 


au. 

&x ; 


yy 


dv. ? 

ay* Y xy " “xy 


au + av 
ay ax 


(43) 


Here b and b are the x and y component of the body force. The 
x y 

material properties E, v, and G have their usual definitions and fig. 4 
defines the coordinate system. The strains may be substituted into the 
compatibility equation to yield 


2 2 2 
o a cr bo n 9 o 

(1-v 2 ) ii - v(l+v) — ^ + (1-v 2 ) 


v(l+v) 


a*o 


ay 

a 2 i 

2( l+v) — 

axay 


ay 


ax 


ax 


XX 

r 


(44) 


The stresses in this equation can be written in terms of a stress func- 
tion F and a scalar body potential V. The body forces b and b are 

x y 


related to V by 


x dx 


h dV 

b y W • 


(45b) 


and are used here, as mentioned In Chapter 2, to represent the pin 
loading from the outer laps. By considering the equilibrium equations, 
the stresses are defined as 


With these definitions, the equilibrium equations are automatical ly 
satisfied. Substitution of these iefinitions into eq. 44 yields 


_4 C l-2v 2,. n 
V F + — 7 V = 0 . 

1 -v 


For all cases where the body forces are derivable from a potential, 
the potential V can be written as 


V ■ V 2 Q . 


(48) 
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Here Q is another scalar function. The redefinition of V further simp- 
lifies the biharmonic equation to 


»V * kQ) = 0 , 


(49) 


where < = — : — . 

1 - v 

It is convenient to take advantage of the biharmonic nature of the 

governing differential equation by using a complex variable approach. 

4 

In this approach the V operator is written in terms of the complex 
variable z and its conjugate 7. Here z and 7 are defined by 


z = x + iy (50a ) 

7 = x - iy. (50b) 


In what follows, the overbar will be used to denote complex conjugate. 
Using eqs. 50, the following differential operators may be defined: 


D, - 


5_ _l ,a_ . d. 

az = 2 l ax ’ 1 ay J 


(51a) 


n d l .a . a . 
D ? = - = 2 ^ax + 1 


az 


ey J 


(51b) 


By using these operators, the biharmonic equation, eq. 49, can be 
written as 


16 D 1 D 2 D 1 D 2 (F ♦ <Q) = 16 — ^ ^ (F + <Q) = V 4 (F + <Q) = 0 (52) 

az az az az 
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Despite the use of complex variables, by eq. 46 F Is real. In 
addition, V and Q are real. The form of F can be found by Integrating 
eq. 52. First the equation Is Integrated once with respect to z and 
then once with respect to z. Integrating with respect to z results In 

~ ~ (F + <Q) = F (7) , (53) 

dz dz dz 

where F^z") Is an arbitrary function of z and any constant of Integra- 
tion Is incorporated Into r^T). Integrating eq. 53 with respect to z 
gives 

2-2- (F ♦ «0) . -F (I) ♦ F (z) . (54) 

dz dz 1 £ 

Here F 2 (z) is an arbitrary function of z, including any constant of 
Integration, and 


'F 




Fj(z) dz . 


(55) 


Realizing 


d__ d_ 

dz dl 



(56) 


is a real operator, and knowing (F + <Q) Is real, the right hand side of 
eq. 54 must be real. Thus v F^(z) and F^(z) must be conjugate functions, 
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i .e 


• • 


F ? (z) = W(z) 

( 67 ) 

F^z) = W(z) . 

(58) 


Then eq. 54 becomes 

— (F + kQ) = W(z) + W(7) . (59) 

sz az 

The equation above can also be obtained by noting that every har- 
monic function (here W(z) Is harmonic) is also biharmonic. In addition, 
the conjugate harmonic function is required to insure that the stresses 
remain real valued. 

Two more integrations are required. Integrating eq. 59 with 
respect to 7 yields 


y (F + <Q) » zW(z) + 'W(z) ♦ F 3 (z ) , 


where F 3 ( 2 ) is an arbitrary function of z and 


'W(z) 



W(z ) dz . 


Finally, Integrating eq. 60 with respect to z leads to 

F + *Q * 7' W ( z ) + z'W(z) + 'F 3 (z) + F 4 (7) , 


(60) 


« • 


(61) 


(62) 
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F^(7) being an arbitrary function of z and 'F-j(z) beihQ defined similar 
to eqs. 55 and 61. Again, making use of the fact that F + icQ is real, 
the right side of eq. 62 must be real. The first two terms are conju- 
gate functions and therefore their sum is real. Obviously "F^ (z ) and 
and F 4 (7) must also be conjugate functions. Equation 62 can now be 
written as the familiar Goursat equation (see ref. 46) for isotropic 
materials. That equation is 


F = j [z<{>(z) + z<t(7) + x (z) + x(7)] - *cQ . (63) 

Here 

♦ (z) = 2 v W(z ) (64a) 

4 >(z) = 2 'W(z) (64b) 

x(z) = 2 'Fj(z) 164c.) 

x (z) = 2 F 4 (z) . (64d ) 


The stress function F has now been expressed in terms of two com- 
plex functions, $ and x> and their conjugates. This, at first, may seem 
to be a step backward, i.e., we have Increased the number of stress 
functions. However it will now be shown that this simplifies the steps 
necessary to satisfy the boundary conditions. 


To Illustrate, consider the case of traction boundary conditions. 

For the pin, a and x are known to be N and T, respectively, on the 
r r re 

boundary. It is convenient to use the fundamental stress combinations 
Introduced by Kolosov (see ref. 46). These combinations are defined as 


4 D l¥ ■ °yy * “x* * 2V 


4 DtD.F = a -a + 2i x 
11 yy xx xy 


Taking the appropriate derivatives produces 


2 2 

+ = 2 ( 4> ' (z ) + 4>' (z )) - 4< -*—^2 + 8 _ 

azaz azaz 


yy xx 


- + 21 T . u = 2 ( z 4'(z) <- x"(0) - -M- 


yy xx 


xy 


az 


These variables are often defined as: 


4>(z) s x' (z) 


$(z) = 0 '(z) 


^(z) = 4»' (z) . 


(65) 

( 66 ) 

(67) 

( 68 ) 

(69a) 

(69b) 

(69c) 


This is merely a change in notation that has been used by several 
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authors. Using eq. 69 leads to following fundamental stress combina- 
tions : 


a + a 
yy xx 


= 2(®(z) + 4>(z)) + (8 - 4 k) ^-9- 


dzdz 


yy 


XX 


+ 2i x vu = 2(z$'(z) + y(z)) - 4k ^-9 . 


xy 


dz 


(70) 

(71) 


The first stress combination is the trace of the stress tensor and Is 
invariant under coordinate transformation. Therefore the stresses In 
polar coordinates remain unchanged, i.e., 

a 00 + °rr = a vv + a xx = + $ (*)) + (8 - 4k) -^-9- . (72) 

DD II yy AX az3z 


The second stress combination is not an invariant quantity and can be 
shown to transform as 


°ee - °rr * 21 Ve = <°yy * °xx + 21 V * ' (73) 

where 

10 n*\ 

a = e . (74) 

It is interesting to note that the magnitude of the left hand side of 
eqs. 68 and 73 are related to the maximum shear stress. 

The stress combinations rewritten in terms of polar coordinates r 


and 9 are 
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°G6 + °rr = 2 ( $ ( z ) + ®( z )) + ( 8 “ 4< ) 


°AA • °rn + 21 X rA = 2cJ 2 (Z$'(z) + 'P(z) - 2< - 2 “^) . 


A- 

azaz 

& 


06 rr 


’r6 


(75) 

(76) 


az 


Subtracting the second equation from the first yields a complex stress 
expression which involves those stresses which are known on the bound- 
ary, namely o rr and x^. The expression is 


o — i x — 4>(z) + $(7) - o 2 (z<j>'(z) + Y(z)) + (4-2k) ^-9- + 2o 2 tc -^-9 
rr re ^ 


azaz 


az 


(77) 


In Chapter 2 eqs. 7 and 21, the stresses on the boundary were 
defined in terms of a complex fourier series, i.e., 

CD 

(o - ix J K . = N - IT = T A. e lkG . ( 21 ) R 

v rr re'boundary k v ' 

Therefore, the right hand side of eq. 77 evaluated on the boundary is 
the complex series involving A k . This equation ties the boundary 
conditions directly to the stress function. 

Having illustrated the convenience of using complex stress func- 
tions and having expressed the traction boundary conditions of interest 
here directly in terms of these stress functions, it remains to find the 
specific forms of 4>, ¥ and Q for the pin. First the form of the body 
force potential, Q, will be found. 
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In Chapter 2 it was assumed that a constant body force acts in the 
x direction. Denote this constant as C. By eqs. 48 and 51, the poten- 
tial V is related to Q by 



dzdz 


(70) 


By definition. 


b 


x 



and 



(79a, b) 


Equation 51a yields 


dV 

dz 


1 (g - 1 £) -7 (C-10) - | c 


ay 


(80) 


Integrating once leads to 


V = \ Cz + g x (z ) , 


where g^(z) is an arbitrary function of z. 
But by eq. 51b, 


dV 

dz 


9,'U) = i C|J + 1 &) =7 (C - 10) 


ay 


(81) 


( 82 ) 


Therefore, Integrating with respect to z gives 
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g^z) = | C z + B, (83) 

B being a constant. Combining eqs. 81 and 83 results in 

V = ^C z + |c z + B . (84) 

Since concern is with the derivatives of V (see eq. 45) the con- 
stant B can be set to zero. Therefore 

V = j C (z + z). (85) 

Substituting this into eq. 78 and integrating once with respect 
to z yields 


4 f - 1 c <“ * 7> f ¥ 2 > ■ 


( 86 ) 


where g 2 (z) is an arbitrary function of z. Integrating once more, with 
respect to z, results in 


2 -2 
1 r,z -z . i 


where 


4Q = j C(|- z + j- z) + v g 2 (z) + g 3 (z) , 
'9 2 (z) - j 9 2 (2) <12 


(87) 


( 88 ) 


and g._(z) is an arbitrary function of z. Substituting for Q in eq. 77 


results in the following expression: 
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°rr ' iT re = + <1>(z) " ° 2 ( z$ '( z ) + *( z ) - z ) + j (2-k) (z+z ) (89) 

Attention can now be focused on the form of the stress functions 
$ and Since the functions $ and ¥ are analytic In the domain of the 
pin, they may each be expanded in Laurent series. These series can be 
written as 


and 


®(z) = l 


n 


V 


*(z) = I b n z n , 


(90a) 

(90b) 


a n and b n being complex constants. If the stresses are to remain 
bounded at the origin of the pin, the negative powers of z must be 
eliminated. This leaves the following form for the stress functions: 


00 00 



(91a) 


(91b) 


ie 

z a re = ro , 



9 = tan _1 (y/x) . 


(92a) 

(92b) 

(92c) 


where 
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These equations define the functional forms of c and v used in the pin 
analysis. The stress combination in eq. 89 can now be written explic- 
itly In terms of z. 

In the following chapter, the radius of the hole in the plate will 
be assumed to be unity. If pin/hole clearance (or interference) is to 
be Included in the solution, it must be accounted for in the analyses of 
the pin. Since the clearance has been defined as +\ in Chapter 2, the 
pin must have a radius 


R = 1 - 


(93) 


An interference fit is obtained by selecting a negative value of \. 

If the functional forms for ¥ and s> from eq. 91 are substituted into 
eq. 89, and the result is evaluated on the pin boundary, z = Ro, the 
result is 


ao 

N - iT = Y A o n 
L n 

.CD 


00 00 



l na R n o n + 
n 


- I 


n = £ 


K D n ~ 2 n 

b n-2 R 0 


(Ro+Ro -1 ) 



(94) 


Here eq. 21 has been used. 

The exact solution for the stress in the pin due to boundary trac- 
tion N - IT = e' Kt can now be determined from eq. 94. This can be done 
by setting A n = 1, n = k, and all other A n = 0 in that equation. To 
solve for the coefficients a n and b n which give the stress functions for 
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this state of stress, it becomes a simple matter of comparing terms of 
like powers of o in eq. 94. Table 1 lists the results of this pro- 
cedure. Using eqs. 91, the coefficients in Table 1 completely define 
the stress functions $ and ¥ anywhere in the pin. Using the fundamental 
stress combinations, the stresses may be calculated anywhere within the 
pin, if so desi red. 

The solution technique discussed In Chapter 2, however, requires 

the displacements as well as the stresses. Specifically, the 

displacements on the boundary are required. For this reason the method 

of determining the displacements will now be discussed. 

To begin the displacement analysis, first recall the constitutive 

relations given in eq. 41. Forming the combination e - e + iy 

yy xx xy 

results in 


yy 


- c 


'XX 


+ Iy 


xy 


1+v r 

= T [ °yy 


- a.. 


xx 


21 T ] 

xy J 


(95) 


The bracketed quantity on the right side of eq. 95 is the second funda- 
mental stress combination, eq. 66. The strains may thus be written In 
terms of simple derivatives of the stress function F, i.e.. 


e 


yy 


xx 



1+v 

E 


[4 



(96) 


In addition, using the strain-displacement relations, the strains can be 
rewritten in terms of spatial derivatives of the displacements, i.e., 
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Table 1 


Coefficients of the Stress Functions <j> and v for the Pin. 


A k 

a n 


b n 


A o - 1 

a o = 

1 

b n - 

0 

A -1 ■ 1 

a l = 

1 

IT 

b n " 

0 

\ * 1 

a l = 

It. 2 . * 
2P. 

b n * 

0 

A_2 3 1 

a 2 3 

j_ 

R 2 

b o " 

-1 

A 2 3 1 

a n = 

0 

b o ■ 

-1 

> 

1 

co 

11 

*— • 

a 3 " 

1 

R 3 

b l ' 

2 

' R 

A 3 * 1 

d n ; 

0 

b l * 

1 

‘ R 

A -k = 1 

a k = 

1_ 

R k 

b k-2 

. lidi 

k_ 2 
r k e. 

+ 

II 

»— • 

a n = 

0 

b k-2 

1 


* _ l-2v 

K — « 

1- V 


> ■ 
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dv . ,Bu dUN 

ay " ax u ax ay J 


l»v 

E 



Using eq. 51, eq. 97 may be written as 


2 a < u - iv > « 

a z az 


Integrating eq. 98 with respect to z yields 


aF 


6(u-1v) " - fy + f> 1 (z) 


( 97 , 


(98) 


(99) 


The function h^(7) in eq. 99 is an arbitrary function of z and Is to be 
determined through the use of the first fundamental stress combina- 
tion. This follows. 

Proceeding in a manner similar to that presented above, the follow- 
ing equation is obtained by using the first fundamental stress combina- 
tion: 


( l-2v)[4 3J-; + 2V] - 2G(t ♦ c ) . (100) 

azaz yy 

Again, the strain-displacement relations and eq. 51 can be used to 
rewrite eq. 100 as 


(l-2v) 


[4 U- t 2V] . 2G [aiu±lvl t . 

azaz az az 


azaz 


(ioi) 
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Making use of eq. 99 and Its conjugate leads to 

2 2 2 

(l-2v) [2 + v] - - + n . ' (z ) - ■~-+h'(z). (102) 

dzdz dzdz 1 dzdz 1 

Combining terms and using the relation between V and Q results In 

2 2 

YU) ♦ V< T ) = 4 ('- v ) 1^= + « • < 103 > 

dzaz azdz 

Using previous definitions, 

Hj'( z) + h^T) = 2(1 —v ) [*(z) + ®(z)] , (104) 

or 

h 1 '(z) = 2 ( 1 - v ) $(z) . (105) 

Integrating eq. 105 with respect to z leads to tne form of h,(F) 
needed in eq. 99, namely 


h 1 (z) = 2(l-v) 4>(z) + D , (106) 

0 being a constant. D is related to a rigid body displacement of the 
pin. Since only the elastic response of the pin is required here, D is 
set to zero. 

Everything is now available to determine the displacements in terms 
of the stress functions. This relation is 
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u - 1 v 


1 

= JS 


[ ( 3-4 v ) <|> 


- z$ 


- * + 


C< 

V 


T 2 


(« ♦ f)3 


(107) 


A simple coordinate transformation yields the displacements In tre 
r,0 coordinate system: 


u 


r 



-1 

a 


(u - iv) . 


(108) 


Witn this result 
within the pin can be 
those of Chapter 4 to 


the tractions and displacements at any point 
determined. These equations will be combined with 
solve the problem. 




Chapter 4 
PLATE ANALYSIS 

This chapter will discuss the development of the solution technique 
for an orthotropic plate. The techniques for analyzing the pin and 
plate are similar, but the material properties fo r each case alter the 
specific forms of the solutions. In addition to this, the plate is not 
simply connected. This complicate': the solution by placing constraints 
on the displacements which did not appear in the pin analysis of Chapter 
3. Fortunately, there is no body force in the analysis of the plate. 
This Is a simplification relative to the Din analysis. The following 
discussion explains the analysis of the plate. 

The plate is assumed to be in a state of generalized plane 
stress. For the plate, the equilibrium equations, the strain-displace- 
ment equations, and the compatibility equations are the same as those 
given for the pin analysis, eqs. 40, 42, and 43. Tne constitutive 
equations for the plate are different than those of the pin and are 
given by 


r 


e 


x 




r o 


r 


°x 1 




c 


y 


► 



o 




0 


0 



(109) 


i 
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where 


E £ plate modulus in the x direction 
x r 

Ey = plate modulus in the y direction 

G = plate shear modulus 

v = plate Poisson ratio, 

xy 


Proceeding in a manner analogous to tnat of tne isotropic case, the 
stresses are written in terms of the stress function F: 


xx 


d 2 F 

^7 


yy 


d 2 F 


( 110 ) 


xy 


d 2 F 

dx5y 


Substituting these stresses into the constitutive equations and defining 
the strains in terms of the stress function, the nontrivial compati- 
bility equation becomes 



4 4 

87 + 1 a F 

dx 2 5y 2 ^y bx^ 


0 . 


(Ill) 


This governing differential equation is referred to as a generalized 
biharmonic equation. The ciiaracteri Stic equation of this differential 
equation is 
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f- / ♦ (J- - 2 ^-) V ♦ f- = 0 . (112) 

x xy x y 

Due to energy considerations, the roots of this equation cannot be real 
valued [48], In addition, the roots will occur in conjugate pairs. 
Therefore, the four roots of eq. 112 can be written as 


p 


1 



^2 



Equation 112 is obtained by assuming 


(113) 


F = F(x + py) • (114) 

If this assumption had been made for the isotropic case and this form 
for F had been substituted into eq. 44, with Q = 0, the result would 
have been repeated roots of p = + i. With this, the argument of F would 
be x + iy or x - iy. lnus the motivation for using complex variables to 
aid in the solution of the biharmonic equation for the isotropic case. 

As in the isotropic case it is convenient to use complex variables 
to aid in the solution of the generalized biharmonic equation. Four 
complex variables can be defined as 


z k = x + p k y, k = 1, 2, 3, 4 . 


(115) 
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By eq. 113 


Z 3 * Z 1 an “ Z 4 = Z 2 


With the definition of z k in eq. 115 ; 


*1 = Yj* + 6 X Z and z 2 = Y 2 z + 6 2 z . 


where 


1-in 1+in. 

y. = 1 and 6 , = k = 1, 2 

■' o ^ o 


Differential operators can be defined as 


D k = ay 


' ^k ax * k = 2 » 3 » 4 • 


The generalized biharmonic equation can thus be written as 


D 1 D 2 D 3 D 4 F = 0 • 


Integrating eq. 121 four times yields the following equation: 


F * 4 (W^x+^y) + W 2 (x+^y) + W 3 (x+^ 3 y) + W^x+^y)) 


(116) 


(117a, b) 


(118,119) 


( 120 ) 


( 121 ) 


( 122 ) 


62 


The stress function F must be real. Thus, by eq. 113, the four 
functions W^, k = 1-4 must be related through complex conjugates. The 
solution for F thus assumes the form 

F - \ («!<*,) * W 2 (z 2 ) * WjtZj) + H 2 (T 2 )) (123) 


This is the Goursat equation for the orthotropic case. 

Equation 123 will play a fundamental role in finding equations for 
the boundary conditions, just as in the isotropic case. It is conven- 
ient at this time however, to introduce some notation which was not 
required for the isotropic case. The following identities will be 
useful when using the Goursat equation: 


a_ 

bx 


bz 


1 bz. 



b _b_ 

by U 1 bz x 



b_ _ i f b_ 
bz " 7 '■bx 



b_ 

sy 



b_ 

bz 


1 

7 


dr'ly) 



- V 

‘ “l* 8^“ 


- b_ . b_ 
^1 bx by 


/ — x b b b 

(“l - l*i) -=- * 1*1 &T - ay 

62 1 

/ - x 6 _ - b b 

(^2 ■ ^ bz^ " ^2 b7 + by 


(124) 
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, - . d d a 

^2 ' ^ - ' u 2 dx " dy * 
oz 2 

Using the definitions of eq. 124, the fundamental stress 
combinations can be found. Recall from Chapter 3 the following stress 
combinations : 


o + o 
XX yy 


o 2 F 

azaT 


yy 


- a + 2i t , = 4 ^—5 
XX xy 


(125a) 


(125b) 


Taking the appropriate derivations of the stress function F, and recom- 
bining eq. 121, the following equations for the stresses are found: 


°xx “i [u l 2 “l" + u z\" + Vv' < 126a ) 

°yy = i [W l" + “2" * V + V' ] (126b) 

T xy = + * ^“l" + ^ 2 “2 ' ‘ 3 (126c) 


where 


u 1 • 
"l 


U * I 

W 2 


a 2 w 1 

^T 

5 2 W 2 

TT 

az 2 


(127a) 


(127b) 
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Unfortunately these equations are not useful in their present form. 

This is because the stress functions and Wp must be found through 

boundary condition requirements. The boundary conditions N-i T act on 
i 0 

the boundary z = e . However, the stress functions are written in 
terms of the variables z 1 and Zp. In other words, the analysis has been 
mapped from the z plane to two planes; z^ and Zp. The boundary condi- 
tions have not been defined in the z^ and Zp planes. To remedy this 
problem, a second mapping is required. This mapping maps z^ to and 
z 2 to Cp. In these mappings, and Cp be forced To transform such 
that at the hole boundary = C 2 = z. Under these transformations, the 
hole in the plate is mapped from an ellipse in the z 1 and Zp planes to a 
circle in the ^ and Cp planes. Figure 7 shows the transformations used 
In this analysis. These transformations require a different approach in 
the development of the traction boundary conditions in terms of the 
stress functions. The relationship between the traction boundary condi- 
tion and the stress function will now be derived. 

First consider a line element, ds, in the z plane. Figure 8a shows 
this element along with its unit normal vector n. From the geometry of 
fig. 8a, the following equations relate dz to the line element ds: 

dx = ds cos (-| + a) 

dy = ds sin (^ + a) 


dz = dx + idy 


.ORIGINAL Pi ■- 
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z plane 
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FIG. 7 Transformations Required for the Orthotropic 
Solution . 
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dz = ds[cos (-| + a) + 1 sin (-| + a)] 

dz = 1e 1a ds 
or 

d 7 = - ie* 1a ds . (128) 

Now consider the tractions T x and T acting on line element ds, 
fig. 8b. These tractions can be related to the normal and tangential 
tractions through the following equation: 

T x - iT y • e" la (N - IT) . (129) 

To find the total force acting on a portion of the circular boundary, 
the tractions acting on line element ds are integrated along the arc 
from the value of s denoting the beginning of the arc to value of s 
denoting the end of the arc. Defining the total forces acting on the 
arc in the x direction as X and in the y direction at Y, the following 
equation is obtained: 


X - iY 



< T * - ’V ds • 


Using eq. 129, the integral car be written as 


(130) 


V 


i Y - / (N - iT)e~ la d$ 


(131) 
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Finally, using equation 128 to write the surface integral in terms 
of 7, the following equation is obtained: 


X 


1Y = i 



(N - i T) dz . 


( 132 ) 


Returning to eq. 124, the tractions can be rewritten in terms of 
the stress function F, by using the foil iwing equation: 


T i * 


ij j 


2 


(133) 


where n, are the direction cosines between the normal to ds and the x 

J 

and y directions. Here 


n 


If and n 2 


ax 

ds ’ 


(134) 


Using eqs. 110 and 128 in eq. 133, the tractions on ds become 

2 2 2 2 

t i t 5 f ay a F ax . r a F ay a f ax \ 

x y = d 2 as axay as axay as ' dx 2 as J * 


(135) 


This equation is obtained when the chain rule is used to differentiate 

^F^ 


T - iT = i — r — - i 
x 'y as '3x 


ay 


i .e. , 





(136) 


"A 
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Using this result in eq. 130, the total forces acting on the arc of 
arbitrary length can be written as: 


X 


iY = 





(137) 


Combining eqs. 132 and 137, a relation between the boundary condition 
and the stress function is obtained, namely, 

2 = f (N - IT) dz . (138) 

Differentiating the stress function with respect to z yields: 

. aw. aw, aw. aw r 

i [6i + 6 ? + Yi — + y? — ] = / ( N - iT ) dz • (139) 

c dz 1 c dz 2 azj az 2 J 

Note that the stress functions and W 2 are functions of z^ and z 2> but 
the boundary tractions are a function of z. Therefore, the stress 
functions must be transformed to the and C 2 planes. This will now be 
done. 

The coordinates z 1 and z ? were previously related to z by 


2 1 * V * V 


z 2 = V * V 


(117a)R 

(117b)R 
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i 0 

On the hole boundary in the plate, z = e = o. Therefore, eqs. 117 
take the following forms on the boundary: 


= YjO + 6j/o 

(140a) 

= y 2 o + 6 2 /o . 

(140b) 


The mappings from the z, to the C, planes and from the z 9 to Co planes 
require ~ z = 0 on tt1e boundary. This motivates the transforma- 

tions from the z^ and z 2 planes to the Cj and C 2 planes to be 


Z 1 = Y 1 C 1 + 6 1 /C 1 


~ >2^2 + ^9^2 


Solving eqs. 141 for and yields: 


/ 2 

+/z. 


Ci - 


V /z l - Vl 

2y, 


/ 2 . 

* 2 V Z 2 ' ^ y 2 6 2 

2 Yo 


(141a) 


(141b) 


(142a) 


(142b) 


71 


The positive signs in front of the radicals have been chosen to insure 
that the mappings put the and C 2 domains outside the unit circle 
(i.e., the mappings give a plate with a hole, not a disk. A disk is 
obtained with a negative sign). 

ac. ac 

The derivatives -r— : - and will be required later. These deriva- 
OZ« o z ^ 

tives are found by differentiating eqs. 142. These derivatives are 


Simi larly , 


azj ( z i = Y i c i + V c i^ 


dz 


z, " c l z l + V 


ac, ac, 

2 y i c i az^ ' z i aiy ‘ c i 


ac 


l 


( 2 y i c i ' z i^ azj = c i * 


5C 2 

(2y 2 C 2 - z 2 ) = C 2 . 


(143a) 

(143b) 


The definitions for C 1 and C 2 which appear in eqs. 142 can be substitu- 
ted into eqs. 143 to yield: 


or 


and similarly 


+ *^i ^ ~ 4y, 6. ac. 


2t. 


ac 

a! 


I 

l 




c 


2 



(144a) 


(144b) 
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Witn these results, it is now possible to define the stress functions of 
eq. 139 in terms of the new coordinates and C 2 » i.e., 


dW. 

(145) 

aw 

8^ E 2 *2< C 2> ’ 


The specific forms of the stress functions ^(Cj) and (t> 2 (C 2 ) will now be 
defined. 

Each component of stress has been defined by eqs. 126. In these 
equations, the second derivative of and are required. Equation 
145 has defined the first derivatives of and W 2 in terms of the new 
stress functions ^(Cj) and 4>,,(C 2 ). Equations 145 must therefore be 
differentiated to find the forms of W^ 1 ' and W 0 ' ' in terms of tne new 
stress functions. This differentiation leads to: 


Simi larly , 


a 2 w 


bz, 


_d 

dz 


dW, 





= 2 


aCj 



dt>2 bC 9 

ac^ az^ * 


(146a) 


(146b) 
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dC dC 

The partial derivatives - — - and have been defined by eqs„ 144. This 

a*. 0Z 1 SZ 2 

leaves and -*-==■ to be defined. Since any analytic function can be 

BC 1 aC 2 

represented by a Laurent series, these derivatives will be represented 
in the following forms: 


5$ ® 

1 r r n 

ac , = 1 a n "“I 



l P n < 



Equations 146 now become 


a 2 w 


az. 


■ 2 n 


JJ 


e 2 w 


l - Vi 


az. 
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(147a) 

(147b) 


(148a) 

(148b) 


The boundary conditions for the plate exist on the inside of the 
hole and at the ends of the plate. The conditions at the inside of the 
hole will be satisfied using the method described in Chapter 2. The 
other conditions are p’aced n the stress distribution at z = ». These 
conditions are: 

1) The stresses at z = ® are bounded, and 

2) Tne plate must remain in equilibrium. 

The geometry of the plate (i.e., infinite dimensions) and condition (2) 
above, require the stresses at infinity to not only be bounded, but they 
must be zero in magnitude. From eqs. 142 it is apparent that 
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if z 1 = ® and = ®, then Cj = « and C 2 = 00 • When infinite values 

of and C 2 are substituted into eq. 148 it is evident that some of the 

coefficients a and B must be zero if the stresses at infinity are to 
n n 

be zero. These coefficients are: 

a = 0 
n 

for n > 0 . (149) 


Pn = 0 


b<f 1 d $ 2 

This leaves the following forms for the derivatives -r-=— and — — - : 

ac 1 ac 2 


a*. 


ac, 



VI 

Mt 2 - l 


1 n=l 


n+L 


(t 2 - l 

’2 n = l 


O-i 
j-i J » 


n 4 -! 


(150a) 

(150b) 


where 


6 o 5 P-1 

- na n 5 r *-(n+l) n > 1 
- nb n 5 P-(n+l) " ' 1 


Integration of eqs. 150 yields the stress functions required in eq. 
139. These stress functions are 


3W 

dz 



2 *i<h> 


2 (‘o' n 


=1 


+ 


CO 


L 

n = l 



(151a) 
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5W °° D 

9zT s 2+ 2 (C 2 ) - 2(b 0 In c 2 + J 77 ) • US'") 

2 n=l Co 

Equation 139 is now fully defined in the C 1 and C 2 planes and takes the 
form 

+ X 2 * 2 CC Z ) *■ Y ! <t> x ( C x ) + Y 2 4> 2 (C 2 ) ] = y (N-IT)dz . (152) 


It is important to realize that the right-hand-side of eq. 152 is inte- 
grated along the boundary arc to some particular value of s. The left 
hand side is evaluated at the value of z = x + iy on the boundary arc 

which corresponds to the end point of the arc of length s. 

i 0 

Recalling that C 1 = C 2 = z = e on the boundary, it is now possible to 

i 0 

define eq. 152 in terms of a common variable a = e . For example, on 
the boundary the stress functions <t^ and become: 


♦i<°> =■ a 0 ln ° * i, -i? 

n = 1 0 

00 b 

4 > 9 (o) = b In o + l — £ 
2 X '0 h 

n = l a 


(153a) 

(153b) 


Differentiating eq. 152 with respect to o (= z on the boundary) now 
yields 


8*. »♦, 8^ ^2 3o 

6 1 6 2 - T 1 do - r 2 80 - " N ‘ T • 

00 oa oa 00 


( 154 ) 


ft 2 

The derivative is equal to -0 , thus eq. 154 becomes 

da 
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3 * 3 * 3 » 3* 2 

6 l 7" * 6 2 7- _ - 0 n aT ' 0 t 2 aT * N • ,T ' 

do do 


(155) 


d$ d<fr 

The derivatives r— and -*4 were defined in eq. 150. Upon substitution 

oQ 2 

of o for Cj and c^, these derivatives become: 


d 4> , , °° a 

J- 1 = 1 (a n - y n -3.) 
do a ( o L . n ; 

n=l o 

d4> « b 

■t— ^ =■ — (b - y n 4) . 

do o v o L . n J 

n=l a 


(156a) 

(156b) 


Substituting eqs. 156 into eq. 155 results in 


- W ‘ n a n°‘ ( "’ 1> ) - *2 (b 0° - I n b n o' (n ‘ 1> ) + 
n=l n= 1 

(157) 

CD CO 

6,(a o - I n a o n + ^) + 6,(b o - ] n b o n + ^) = N - IT . 

1' o i n ' 2' o „**. n ' 

n=l n=l 

Just as in the isotropic case, N - i T is expanded in a Fourier 

series and like powers of o are compared to obtain the coefficients a n 

and b . This procedure is summarized in Table 2. The coefficients a , 
n o 

b Q , a Q , and b Q require special consideration. This is as follows. 

Comparing the coefficients multiplying o + * yields two equations for 

the four unknowns a„, b rt , ? and ¥ . To obtain two more equations, a 

o o o o 

restriction on the displacement field in the plate must be imposed. 
Mathematically, the logarithmic terms in eqs. 151 are multivalued. If 
this were left as such, the displacements in the plate would also be 
multivalued. It is necessary, therefore, to impose restrictions on the 
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Table 2 

Coefficients ror tbe Stress Functions $ and <j> 2 . 


A n 


a n 


b n 

A o 

a l = 

5 1 

^^1 " 6 ? y 2 ) 

»1 = 

f>2 

l ^ y 2 " V 

A -i 

a 2 = 

6 2 

1/C2(T 2 -sJ-Vj)] 

bo = 

im 

l/C2( Tl -^T 2 )] 

A .i 


★ 


★ 

A -2 

a 3 = 

» 2 

wi'2 - sf n> 3 

b 3 

l/[3( n -^r 2 )] 

a +2 

a l = 

i/(6 2 ^- 6i) 

b l * 

l/(6l^- 6 ? ) 


A -k a (k +1) * '/[( k+1 )(T 1 '5jv 2 )] b (k*l) ’ 1 /[( k+1 )(''i'6^2> ] 

A *k * (k -l) ' 1/C(k-1)(6 2 -^ -8j)] b (k l) - l/[( k -l)(6^ -6 2 )] 

*See text for explanation 


coefficients a Q and b Q which will force the displacements to be single- 
valued. In this process, two more equations involving these coeffi- 
cients are generated. There are then four equations with which to solve 
for a Q , b 0 , a Q , and b Q . The equations needed to enforce the single 
valued displacements will now be developed. 

To begin, it is necessary to develop an equation which relates the 
displacements to the complex stress functions. The procedure for the 
orthotropic case follows that of the isotropic case exactly. Specific 
differences arise due to the forms of the constitutive relations for 
each case. Due to the complexity of the algebra involved in the ortho- 
tropic case, the derivation of the displacement equations will not be 
presented here. Instead, an outline of the method is presented below. 

1) Write the strains in terms of the displacements 

5u 8v au av 

e x "ax e y ay Y xy "ay ax * 

2) Use the stresses of eqs. 126 in the constitutive relations to 
find the derivatives of the displacements in terms of the stress 
functions W^Zj) and W 2 (z 2 ). 

3) Combine these equations and integrate with respect to z to 
obtain the displacement equation. 

The following equation is the result of this procedure: 

4 ( u + iv) = p^' + p 2 W 2 * + q^' + q^' , (158) 


where 
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P 1 E 


. + 1- u 2 + 1- f- 1 - _ u 2 ) 

E E P 1 p, [ E E ^1 > 
xx ^1 y x 


P 2 ^ 


V xy 1 2 1 r l jg_ 2^ 

e + r ^2 + 77 % " E ^2 ) 

xx 2 y x 


9i = 


99 E 


xy .l-2i r l _xy_ - 2^ 
■ E + T P 1 + =~ 1“ ' F Pi J 
x x 


l l 


xy 1 - 2 i r _l_ _x£ - 2 >, 
"E + E ^2 + - ‘ f P? J 

x x 


Recall from eqs. 151 the following definitions: 


aw. 


cd a 


V ■jr = ! [‘."'i t n I, TTf) 

1 n*l C 1 

dW °° b 

W 2* = 3zT = 2lb o ln C 2 + E ~n^ ' 
2 n=l C 1 


(159a) 

(159b) 

(\59c) 

(159d) 

(160a) 

(160b) 


i 0 

When transversing the boundary z = e from 0 = 0 to 0 = 2n, it is 
evident that each term of the stress functions W^' and W^' returns to 
their original value except the logarithmic terms. These terms increase 
by a factor of 2ui. The displacements at 0 = 2n and 0=0 are therefore 
related through the following equation: 


(u 4 iv > 9 = 2 „ = 2 ,li leiV p 2 b o- q iV q 2 F o> 4 (“ 4 iv >e=o (16l) 


Thus, if the displacements are to be single valued. 


and its conjugate 


Vo * p 2 b o - Vo - Vo * 0 • (162b > 

Using these two equations and the equations obtained by comparing like 

powers of a in eq. 157, a , b , 7 and b can be found. This will be 

oooo 

done numerically in the computer program developed for this study. 

With the analysis of the orthotropic plate complete, the results 
can be combined with the reults of Chapter 3 to implement the solution 
scheme described in Chapter 2. 


Chapter 5 

. .TAILS OF THE NUMERICAL SCHEME 

As described in Chapter 2, the various coefficients multiplying the 

unknowns A^ in eq. 39 are computed by the closed-form complex variable 

techniques described in the previous two chapters. The constants a^ and 

t> k tabulated in Tables 1 and 2 are used to compute the stresses and 

i k O 

displacements of the plate and pin due to the traction e . Once all 
the quantities have been computed, they are assembled into the system of 
linear equations, eq. 39. The system is solved using a standard linear 
equation solver. 

Modifying the system of equations 

Operationally, eq. 39 had to be modified in order to effect a 

solution. As they are stated, the system of equations in eq. 39 cannot 

be solved. The reason is as follows: The elasticity solution for the 

displacement of the pin and the displacement of the plate to the trac- 
i k 0 

tion N-i T = e automatically yields 

(u ) k « (u ) =0 <a 9 = 0 , k = 0, +1, ... + N . (163) 

pin plate 

This is due to the antisymmetric nature of the functional dependence 

on 0 of the tangential pin and plate displacements to traction N-1T 
i k 0 

= e . Also due to the nature of the solution, the conditions 

(T r0 )k = ( T r0^k = 0 at 0 = 7i, k = 0, +1 + N . (164) 

pin plate 
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are automatically satisfied. Tne collocation procedure included the 
points 0 = 0 and 9 = it. This is necessary because there is one other 
condition at each point that must be enforced. These conditions are 

(u ) . = (u ) . @0=0 and o =0 @ 0 = n. (165, 166) 

v r'pin ' r' plate rr v ' 

The conditions 

u Q - u Q =0 @0=0 and (t ) =0 @ 9 = n( 167 , 168 ) 

u • y i . ry 

pm plate 

were eliminated from eq. 39. To include them would result in two rows 
being identically zero and producing a singular matrix. Thus the number 
of equations solved and the number of A k 's was 2N-2. 

Finding a and p 

As has been mentioned previously, actually solving the reduced set 
of linear equations required having numerical values for a and p. A 
specific set of values were chosen for a, p, and 6, the rigid body pin 
displacement. The system of equations was then solved. In all likeli- 
hood the values of a and p chosen were not the correct values for the 
specified value of 6. As a consequence, one or more of the physical 
realities of the problem, expressed as what will be termed auxiliary 
boundary conditions, were violated. Subsequently, an iteration process 
was used to find values of a and p which satisfied all the conditions. 
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The iteration processes was divided into two parts, 1) the finding 
of p, and, 2), the finding of a. The two parameters were not found to 
be strongly dependent on one another. 

The key to the Iteration process for finding a and p was the satis- 
faction of the auxiliary boundary conditions. These auxiliary boundary 
conditions are the combination of a restatement of the conditions stipu- 
lated by eqs. 1-6, and statements of the physics of the problem. 
Specifically, the auxiliary boundary conditions are: 


N = 

T = o @ e = p 

(169) 

N = 

T = 0, P < 9 < n 

(170) 

N < 

0 all 0 

(171) 

T = 

pN @ 9 = a 

(172) 

T < 

pN 0 < 9 < a 

(173) 


Equations 169 and 172 are written explicitly because these were the most 
powerful Indicators of convergence to the correct values of a and p. 
These above five equations express the fact that the normal traction 
must remain compressive everywhere and that it must be zero at the value 
of 0 corresponding to the loss of pin/plate contact. Also, in the 
region of no slip, the friction-induced shear traction, must be less 
than or equal to the normal traction. At the end of the no-slip region 
the shear traction must be exactly equal to p times the normal traction. 

Figure 9 diagrams the procedure used ir. the iteration process to 
determine the half-contact angle, p. In this figure, a typical stress 


- 1.0 


3 Too Small 



Fig. 9 Iteration Procedure for the Half-Contact Angle, b, 
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distribution for tne plate at R = 1 is shown for 0° < 0 < 90°. Tne 

normal stress, o , in the region of 0 = 0 is shown for three choices 
rr 

of 0. As shown In fig. 9a, if 0 is too great, the normal stress becomes 
tensile -- a physical impossibility in che contact region. If 0 is too 
small, however, as shown in fig. 9c, the normal stress is compressive 
but not zero as required by the boundary condition at e = 0. The cor- 
rect choice of 0, fig. 9b, yields a zero value for the normal stress 
at 0 = 0. All cases considered nave followed this description of the 
iterative procedure exactly, and each solution appears to be unique. 

A similar iterative procedure for determining the region of no- 
slip, - a < 6 < a, is shown in fig. 10. In this case, the shear 

stress, x , at 0 = a is analyzed. If the shear stress at e = a is 
re 

greater than |nc rf .|, as in fig. 10a, then the no slip angle is too 

large. On the other hand, if < jno rr |, as in fig. 10c, then the no , 

slip angle is too small. Wnen the correct choice of <i is made, as In 
fig. 10b, then the boundary condition at 0 = a (i.e., T pe * |po rf .|) is 
satisfied. As in the case of the contact angle, this procedure was 

t • 

found to be general, and all solutions to date appear to be unique. 

The procedure for finding a and 0 was to set a to a value 
corresponding to the angular location of the first collocation point 
away from 0=0. With this value of a, the value of 0 was found by the 
iteration method discussed in fig. 9. Once a value of 0 was found, the 
Iteration procedure focused on finding the va’ue of a. When this was 
done, 0 was checked. There generally was little coupling between the 



two variables. 


c*i XD^ot-j 



Fig, 10 Iteration Procedure for the No-Slip Angle 
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For lower values of the coefficient of friction (e.g., n = 0.2), 
the no-slip region was found to be small. Often times the angle a was 
less than the spacing between 0=0° and the first collocation point. 

In these cases, the no-slip boundary condition was satisfied at only one 
point i.e., 0 = 0°. Larger values of the coefficient of friction 
(e.g,, \i = 0.4) enlarged the no slip region. For these cases, the 
number of points at which the no-slip boundary conditions was satisfied 
were increased. Due to the limited size of the no-slip region, it was 
found to be very simple to iterate to the correct value of a. The half- 
contact angle, p, on the other hand, varied greatly from case to case 
and generally took more iterations to converge to the correct solution. 

Collocation po i nts 

Overall, however, the convergence of the solution was found <;o be 
dependent on the location and number of points used to satisfy the 
boundary conditions (i.e., the collocation points). As mentioned in 
Chapter 2, two boundary conditions were enforced at each point. In 
earlier stages of programming, each boundary condition was enforced at a 
different, point. This arrangement proved to have poor convergence 
characteristics and was abandoned in favor of enforcing two boundary 
conditions at each point. 

Additional convergence problems were encountered when the spacing 
between the collocation points was not uniform. The analysis requires 
that the number of points in the contact region and those outside the 
contact region be specified separately. It is quite possible that the 


C 
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point spacing In one region could be quite different than in the other 
region. In such cases, the solution was found to oscillate in the 
region with the widest spacing. This problem was easily overcome, 
however, by ensuring that the point spacing was as uniform as possible 
for the entire region 0 < 6 < n. 

It was found that if fewer than ten points were used to satisfy the 
boundary conditions, the solutions were inaccurate. Twenty to thirty 
points along the boundary gave reasonable results. Forty points were 
used to obtain the results presented in this report. Forty points was 
close to the upper limit of the interactive capability of the computer 
used in this analysis. 

The solutions presented here all have 78 unknown A^'s in the series 
representation of the boundary traction. The 78 terms were divided 
evenly between negative k terms and positive k terms. The solutions 
were forced to satisfy a total of 78 boundary conditions at 40 points 
around the half-circle. This means the boundary conditions are met at 
roughly every 4.6° around the circumference of the hole. 

Computer programs 

The computer programs were written at NASA Langley Research Center, 
Hampton, Virginia. The machines that were available for this work were 
Control Data Corporation Cyber 1 7 3 ' s . Control Data's Fortran, version 
5, language was used throughout the analysis. Due to the finite memory 
limitations of these machines and the extensive use of complex variables 
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in the analysis, it was necessary to use single precision accuracy. 
Single precision accuracy for the Cyber 173 is 14 digits. 

Run times were typically 12-15 seconds. This made interactive 
analysis practical. In fact, one of the goals of the numerical scheme 
was to make it efficient enough so that convergence to a and p could be 
effected by means of interactive computing. A "human In the loop" made 
the scheme more attractive than imbedding the stress analysis in an 
iteration loop with convergence criteria. With graphics and short run 
times, the program was ideal for interactive usage. With some practice, 
a solution could be found in as little as four iterations. 


Cnapter 6 
NUMERICAL RESULTS 

Comparison with past work 

To build confidence in the present technique, numerical results 
obtained from it were compared with the results of other researchers. 

The results of previous researchers generally include more restrictive 
assumptions thin were used in this study and were often more restricted 
In overall scope. However, they from a basis with which to compare the 
current results. 

The studies which are most similar to the present work were pro- 
duced by de Jong [3, 11] and Crews [16]. Recall that de Jong's work 
represented a rigid pin in an infinite orthotropic plate. This work 
Included friction and only considered loading, not unloading. In addi- 
tion, de Jong studied only the linear, push-fit case. Crews, on the 
other hand, studied elastic pin-loaded plates of various widths. His 
study did not include friction and the clearance was not varied. None- 
the-less, the papers by de Jong and Crews offer good results for compar- 
ison with the present study. One of the comparisons *s presented in 
figs. 11-12. This comparison is for the frictionless case of a rigid 
pin-loaded isotropic plate, Ine contact stresses and the hoop stress 
are shown for the region 0 < 0 < 90° in these plots. Symmetry and 
antisymmetry conditions can be used to find the stresses in the region 
-90° < 0 < 0. Outside these regions the contact stresses are zero and 
the hoop stress decays stecdily to a near zero value at 0 = 180°. This 
display of results will be used throughout the rest of this report. 
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Figure 11 shows a comparison of the normal stress, o rr , as computed 
by de Jong [3], Crews [16] and the present investigation. Figure 12 
shows a comparison of the hoop stress, a . As can be seen, there Is 

v7 U 

excellent agreement among the three studies. There is a difference 
among the three, however, at the end of the contact region (9 = p). 

This is not surprising since neither de Jong nor Crews Iterated to find 
the value of p. In contrast, iteration to find the correct value 
of 0 was a key step in the present analysis. 

Another comparison with previous results is shown in figs. 13 and 
14. These figures display the contact and hoop stresses as a function 
of circumferential location. The results are for a rigid pin-loaded 
quasi -Isotropiu plate. Friction is included in both de Jong's [11] and 
the present study. The agreement between the two solutions Is good 
except in the no-slip region, near (9 = 0), and at the end of the con- 
tact region. Recall, de Jong assumed boundary conditions In the no-slip 
region. Also, de Jong did not iterate to find p. The boundary condi- 
tions were not approximated in the present study, hence differences 
might be expected. The following discussion demonstrates the analysis' 
flexibility by presenting the results of a limited parametric study. 

This parametric study was designed to show the influence of the various 
parameters on the response of the plate. As such, the parametric study 
presents new and interesting results regarding the stress In pin-loaded 
orthotropic plates. 
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90.0 


Fig. 12 Comparison of the Hoop Stresses for an Isotropic 
Plate Loaded by a Rigid, Frictionless Pin. 

(Crews - ref. 16, de Jong - ref. 3) 
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Fig. 14 Comparison of the Hoop Stresses for a Quasi - 
Isotropic Platf Loaded ey a Rigid Pin with a 
Friction Coefficient of 0.2 (de Jong - ref. 11) 
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Parametric stud ies 

A standard ca r e has been adopted to show which parameters have the 
greatest effect on the problem. The parameters for the standard case 
are 


pin - steel , 

friction coefficient - 0.2, 

loading direction - positive (i.e., ^ = -0.2) 

clearance - \ =0.01, 

rigid body pin displacement - 6 = 0.07. 

A steel pin (E = 30 x 10 6 psi , v = .25) was chosen to simulate steel 
bolts commonly used to join composite panels. Various values of the 
coefficient of friction between steel and graphite epoxy have been 
found, but u = 0.2 seems to be the most widely accepted value and will 
therefore be used as the standard case in this report. A clearance of 
1% has also been found to be reasonable for typical joints tested in a 
laboratory environment. The parametric study will begin by varying the 
plate properties. Various laminates of graphite epoxy will be used to 
study the effects of plate properties on the stress distributions. The 
properties were chosen from the report by Crews [16], They are listed 
in Table 3. As can be seen, these laminates represent a wide variation 
in plate properties. 

Tnese properties were used along with the other parameters in the 
standard case to compute stresses. The numerical results appear in 
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Table 3 

Plate Properties 


Laminate 

E x • 10 " 6 (psi) 

E y • 10 " 6 (psi) 

v G y w • 10‘ 6 (psi ) 

xy xy \r i 

[0/+45/90] s 

8.40 

8.40 

0.310 

3.20 

[0] 8 

21.31 

1.58 

0.380 

.93 

C90] 8 

1.58 

21.31 

0.028 

.93 

[0/90J 2s 

11.53 

11.53 

0.052 

.93 

C±45] 2s 

3.26 

3.26 

0.735 

5.47 

[0/+45] s 

9.31 

4.09 

0.699 

3.96 
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figures which occur In pairs according to the plate material. The first 

figure of each pair gives the normal and shear contact stresses, 

a and and the hoop stresses, o , for the plate at the hole 
rr rw 04 

edge. This Is similar to the presentation of results for comparison 
with other Investigators. The second figure of each pair gives the net- 
section and shear-out stresses for each plate. The net-section and 
shear-out planes are shown in fig. 15. The stresses along these planes 
are typically studied when the prediction of the failure mode is 
required. Often In design studies, the net-section failure mode is 
suppressed. This, however, can lead to shear-out failures. The 
stresses along each plane are therefore presented in an attempt to 
describe the overall response of each layup. In each figure the 
stresses have been nondimensl onal ized by F, which for a plate of unit 
thickness. Is the bearing stress on the hole. In addition, since the 
hole In the plate has a unit radius, i.e., equals 1.0 in., the parame- 
ters 6 and \ are effectively nondimensional i zed by the hole radius. 

Figure 16 shows the stresses at the hole-edge for a quasi -1 sot r op1c 
laminate. Th<s laminate is quite popular In current applications of 
composite materials. As can be seen, the hoop stress is compressive 
at 6 = 0 and becomes tensile as 0 increases. The maximum hoop stress 
occurs near the end of the contact region. The negative hoop stress 
at 0 = 0 Is due to the effects of friction. The net-section and shear- 
out stresses, shown In fig. 17, are as could be expected. It is inter- 
esting to note that the shear-out stress peaks within a hole radius of 


the net-section. 


99 



Fig. 15 Net-Section and Shear-Out Planes 
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y/R 

a) Net- Section Stress 



x/R 

b) Shear-Out Stress 


Fig. 17 Net-Section and Shear-out Stresses for a [ 0 / ±£ 45 / 901 s 
Laminate Under Standard Case Conditions. 
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Figure 18 shows the contact and hoop stresses for a [0] Q laminate 
while fig. 19 shows the net-section and shear-out stresses. This case 
Is an extreme example of the effects of orthotroplcity. There are 
several unusual features for the stress distributions illustrated 
here. One of the most prominent features is the distribution of the 
normal stress, o rr . Many analyses assume that a sinusoidal distribution 
is valid independent of laminate. It is clear from fig. 18 that such an 
assumption for a highly orthotropic laminate is in error. Another 
prominent feature is the peak in the hoop stress near 0 - 90°. Not only 
is the magnitude of the peak greater than that of the quasi -i sotropic 
case, but the peak occurs outside the contact region. It is interesting 
to note that the hoop stress changes sign several times within the 
contact region. The net-section stress for the [0]g laminate, fig. 19a, 
looks similar to the net -section stress of the quasi-isotropic laminate, 
but the magnitude at the hole edge is considerably greater. On the 
other hand, the spatial dependence of the shear-out stress is consider- 
ably different than the quasi -isotropic case. The magnitude of the 
shear-out stress for the [0] 8 laminate is less, but the shear strength 
of this laminate is considerably less than that of the quasi -isotropic 
laminate. It can be expected that a [0]g laminate would fail in shear- 
out. 

Figures 20 and 21 document the stresses for a [90 ]q laminate. This 
laminate represents the other extreme in orthotropicity. In practice, 
such a laminate would never be used. The noop stress in fig. 20 has two 
peaks of roughly the same magnitude. The hoop stress shows very little 
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y/R 

a) Net- Section Stress 



x/R 

b) Shear-Out Stress 


Fig. 19 Net-Section and Shear-Out Stresses for a [o] g 
Laminate Undf.r Standard Case Conditions. 





106 



2 3 4 5 


y/R 

a) Net- Section Stress 


.75 



b) Shear-Out Stress 

Fig, 21 Net-Section and Shear-Out Stresses for a [90] g 
Laminate Under Standard Case Conditions. 
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similarity with the previous two cases. Note that o Q „ is compressive 

over about half the contact region. The normal stress, a , is 

rr 

radically different than the previous two cases. A sinusoidal approxi- 
mation would be very much in error for this case. The shear-out and 
net-section stresses for the [90]g laminate are nominally the same as 
the quasi -isotropic case (fig. 17). However, the net-section stress Is 
no doubt tne failure producing stress for the [90]g laminate. 

Figures 22 and 23 illustrate the stresses for a [0/90] 2 S lami- 
nate. This laminate, though not as popular as the quasi -1 sotroplc 
laminate, does have some applications. As with the [90]g laminate, 
there are two peaks in the hoop stress. The magnitude of the peak 
nearest 9=0, however, is less than the peak near 9 = 90°. This is no 
doubt due to the influence of the 0° laminae. As with the [0]g case, 
the maximum hoop stress occurs outside the contact region. It siiould be 
pointed out that the maximum hoop stress at 9 = 90° is not too much less 
than the maximum stress for the [0]g laminate. This is due to the fact 
that the 0° laminae, rather than the 90° laminae are responsible for 
transmitting the load past the net-section. Thus, similarity with the 
[0]g laminate might be expected in this region. The net-section and 
shear-out stresses, for the [0/90] 2s laminate, fig. 23, show ,io 
unexpected results. 

Figures 24 and 25 show the numerical results for a [±45]^ s lami- 
nate. Such a laminate can be used where high shear strength and/or 
stiffness is required. Tailoring studies [34] have suggested that, for 
higher performance joints, the region near the pin should contain a high 
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v/R 

a) Net- Section Stress 



Fig. 25 Net-Section and Shear-Out Stresses for a [-^5] 2S 
Laminate Under Standard Case Conditions. 
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percentage of i45 ply s . The most obvious feature of the hole edge 
stresses In the [+45^ laminate Is the maximum hoop stress occurring 
at 6 ■ 45°. This should be put In context because the maximum strength 
of the laminate is also at 0 = 45 w . The magnitude of the peak stress 
at 0 a 45° is comparable with the peak stress at the net-section of the 
quasi -1 sotropic laminate. The other components of stress for the 
[+45]2 S laminate are similar to that of the quasi -isotropic case. 

Finally, figs. 26 and V show the stress distribution for a 
[0/+46] $ laminate. This laminate has many applications. The differ- 
ences between the results for this laminate and the quasi -1 sotropic case 
are small. The hoop stress, however, is higher near 0 = 45° for the 
[0/+45] $ laminate than it is for the [0/+45/90] s laminate. This may be 
due to slightly higher percentage of 45° laminae in the [0/+45] s case. 

The results presented in these figures have shown that the plate 
properties have a significant effect on pin-loaded joints. Other para- 
meters affect the stress distributions as well. Therefore, the plate 
properf.es will now be held constant and the pin modulus, friction 
coefficient, and clearance will be varied. A [02/+45] s laminate was 
chosen for this portion of the parameter study because it represents a 
moderately orthotropic material which is commonly used in the aerospace 
Industry. The properties for this layup were selected from de Jong's 
report [11]. This was done because de Jong presented results of a 
parametric study which complement the current research effort. The 
properties used by de Jong were 
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Fig, 26 S' 
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r.ffp 



y/R 

a) Net- Section Stress 



Fig. 27 Net-Section and Shear-Out Stresses for a [0/ ±45 ] s 
Laminate Under Standard Case Conditions, 


115 


E = 11.61 x 10 6 psi 

E y = 3.43 x 10 6 psi 

v = .662 (174) 

xy ' ’ 

G = 2.40 x 10 6 psi 
xy 

The format used to illustrate the effects of pin modulus, clear- 
ance, and friction will be the same as has been used in presenting all 
results so far. Here, however, each situation studied will involve four 
figures. This was done primarily to emphasize the effect of these 
variables on certain stresses and to avoid cluttering the figures. The 
first two figures of each set illustrate the effect of parameter varia- 
tion on the norma! and the shear stresses. The second figure illu- 
strates the effect on the hoop stress. The third and fourth figures 
show the effects on the shear-out and net-section stresses. 

Using the laminate elastic properties given by eq. 174 along with 
the values given for the parameters of the standard case, the pin modu- 
lus is now varied to see what effect this has on the stresses. Three 
values of pin modulus are studied: a perfectly rigid pin, a steel pin, 
and an aluminum pin. Tne stresses tor these cases are presented in 
figs. 28-31. It is evident from these plots that the pin modulus plays 
a very small role in this problem. Apparently a rigid pin assumption 
can be made without introducing any significant error in the problem. 

In figs. 32-35, the standard case with variable clearance, \, is 
shown. In these plots, the rigid body pin displacement is held cor - 
The nonlinear nature of the problem under these conditions 1u 
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0 22.5 45.0 67.5 90.0 

9 (degrees) 


Fig, 28 Contact Stresses at R = 1 for the Standard Case 
with Variable Pin Moduli i. 





Rigid Pin 

Steel Pin 

Al. Pin 


67.5 9 

6 (degrees) 

Fig, 29 Hoop Stresses at R = 1 for the Standard Case With 
Variable Pin Moduli i . 
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easily observed. In fig. 32 it is obvious that the contact angle 
decreases with Increasing clearance. This causes the location of the 
maximum hoop stress in fig. 33 to vary for these cases. This may have a 
direct effect on the Initial damage location in composite panels. 

The coefficient of friction, n, is varied in figs. 36-39. The sign 
of the coefficient is negative, indicating monotonic leading for each 
case. From these figures It is apparent that the magnitude of the 
coefficient of friction plays an important role in the stress distribu- 
tion of pin-loaded joints. For example, the hoop stress at tne origin 
actually changes size due to the presence of friction. Also note, the 
contact angle, 0, changes very little, however, the no-slip angle, 
a, changes substantial ly for these cases. 

A comparison of interference, push, and clearance fits is shown in 
figs. 40-43. As witn the case described in figs. 32-35, the value of 
the clearance parameter causes the maximum hoop stress to change loca- 
tions on the inside of the hole. 

Finally, in figs. 44-47, the load direction is varied. In these 
figures, n = +0.2 corresponds to a joint which has been loaded past 6 = 
0.07 and monotoni cal ly returned to this value of pin displacement. The 
negative value of n indicates a panel that has been monotoni cal ly loaded 
from 6 = 0 to 6 = 0.07. Load reversal is seen to have a significant 
effect on the stress distribution in the plate. Obviously the noncon- 
servative nature of friction tends to lock in the stresses associated 
with the higher load level. As has been obvious, the effects of clear- 
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Fig. 40 Contact Stresses at R = 1 for the Standard Case 
with Interference, Push, and Clearance Fits. 
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Fig. 44 Contact Stresses at R = 1 
with Unloading, Zero Frict 
Conditions. 
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0 22.5 45.0 67.5 90.0 

0 (degrees) 


Fig. 45 Hoop Stresses at R = 1 for the Standard Case with 
Unloading, Zero Friction, and Loading Conditions. 










ance, friction, and pin elasticity on the net section stress and the 
shear-out stress, for all cases considered is negligible. 

The cases presented in this chapter are a sample of the analysis's 
capability. These figures have demonstrated the flexibility of the 
analysis and give a limited study of the parameters which effect the 
stress distribution in pin-loaded joints. This study should help define 
the importance of these parameters when future design analyses are made. 

Isotropic plate 

The Appendix presents the analysis and numerical results for the 
case of an isotropic rather than an orthotropic plate. Although not 
directly related to the analysis of composites, the pin-loaded isotropic 
plate is an interesting mechanics problem. Hence the results for this 
case, which were derived during the course of this investigation, are 
presented. In fact, the isotropic plate can be used, with a high degree 
of accuracy, to represent a quasi -i sotropi c plate. 




Chapter 7 

CONCLUSIONS & RECOMMENDATIONS 

The work presented In this thesis has been directed towards a 
better understanding of the stresses in composite pin-loaded joints. 

The results of some simple cases were compared with those of previous 
researchers. More Importantly, new results which represent an advance 
In the state of the art were presented. From this analysis, the 
following conclusions are drawn. 

1) The algorithm developed has been shown to be accurate and 
efficient. 

2) The material properties of the plate in pin-loaded joints 
strongly influence the stress distributions. 

3) The magnitude of the coefficient of friction has a large effect 
on the stress distribution in the plots. 

4) Loading direction greatly alters the stresses. 

5) Clearance (or interference) can cause the location of the 
maximum hoop stress to vary. 

6) Pin elasticity has almost no effect on the stress distributions 
for the cases studied here. 

Altnough this analysis has produced a better understanding of the 
state of stress in pin-loaded joints, more work can be done in this 
area. Specifically the following recommendations are made: 

1) Study finite width and end effects. 

2) Study multiple-pin configurations. 
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Develop a failure theory to use In conjunction with the stress 
analysis. 

Include the cases where the loading direction and the prlnclpl 
material directions do not coincide. 

Possibly reformulate the problem so that the pin load Is the 
Independent parameter rather .han the pin displacement. 

Verify the results with an experimental investigation. 
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Appendix 

ISOTROPIC PLATE ANALYSIS 

The research presented in this thesis is primarily concerned with 
the development of an algorithm which predicts the stress distribution 
in pin-loaded orthotropic plates. The mathematics change however when 
the plate material is isotropic. Quasi -i sotropic laminates can be 
treated successfully with an orthotropic analysis, but to remain rigo- 
rous in the mathematical development, the truly isotropic case must be 
treated separately. The mathematics for this case is presented in this 
appendix to complete the analysis of pin-loaded joints. 

The case of an elastically pin-loaded isotropic plate was actually 
developed before that of the orthotropic plate in the current research 
program. This was done because tne analysis of the isotropic plate 
parallels that of the isotropic pin. None of the complicated mappings 
required for the orthotropic case are required for the isotropic case. 
Witn the simpler analysis of the isotropic case, the problems associated 
with the analysis method described in Chapter 2 were solved before the 
orthotropic case was undertaken. Enough differences exist between these 
cases, however, to warrant a separate description of the isotropic case 
in this appendix. 

The pin analysis technique remains the same for pin-loaded isotro- 
pic plates. The same is true for the overall method of analysis 
discussed in Chapter 2. The only difference in the problem formulation 
is the analysis of the plate. Thus, this appendix will be concerned 
only with the development of the analysis for an isotropic plate. The 
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governinq equations for the pin and plate are similar, therefore an 
abbreviated version of their development will be presented here. 

The constitutive, compatibility, and strain-displacement equations 
remain the same as they appear in Chapter 3. The equilibrium equations, 
however, contain no body force for the analysis of the plate. This 
similifies the governing biharmonic equation to the following form 



a 4 r 

dx 2 dy 2 



0 


(Al) 


Integrating eq. Al yields the following Goursat equation 


F 


* \ [z$ + Z4 


x] 


(A2) 


Note that eqs. 47 and 63 of Chapter 3 become eqs. Al and A2 when the 
body forces are set to zero. 

The fundamental stress combinations can now be found in terms of 
the two stress functions. T h ese combinations are defined by eqs. 65 and 
66 of Chapter 3. With the body forces set to zero these equations 
become 


V 


°yy + °xx = + 


(A3) 


a -a + 2ix = 2[z$' (z) + ¥(z)] , 
yy xx xy 


(A4) 


* 
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where 4^(2 ) s x' (z ) 

*(z) = *'(*) 

4>(z) = 4> ' (z ) • 

Using the appropriate transformations, these combinations can be written 
In the r-9 coordinate system as 

°96 + °rr = 2 ^ z ) + 
o 00 - o pr + 2iT r0 = 2o 2 [l4>'(z) + Y(z)] , 

. i9 

where a = e . 

Combining eqs. A5 and A6 yields the required traction boundary condi- 
tion: 

N-iT = $(z) + $(T) - o 2 [7$(z) + ur(z)] . (A7) 

The form of the stress functions must now be defined. As before, 
these functions are expanded in Laurent series 

CD 

»(z) = l ( A8 ) 

— CD 
CD 

’?(Z) = l b n z n . 


(A5) 

(A6) 


(A9) 
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In Chapter 4, the infinite orthotropic plate was required to satisfy the 
following conditions: 

1) The stresses at infinity remain bounded, and, 

2) The plate must remain in equilibrium. 

When these conditions are enforced for the isotropic plate, the follow- 
ing restrictions are imposed on the form of the stress functions 


a n = 0 n > 0 


(A10) 


b = 0 n > 0 
n 


(All) 


Thus the following forms for the stress functions remain 


00 (X 


*(t) - I -a 

n = l z 


*<*> - i -J. 

n=l z 


(412) 

(A13) 


Substituting these forms into eq. A7 and noting that z = a on the bound- 
ary, the following equation is obtained 


oo (X 


N-1T = y A c n = 2 a o' 1 + I T o" + n 2 -2L - b.a - 2 -V . 

n n=l n n=l " n=l a" 1 n^O o" 


(A14) 


As before, like powers of o are compared to find the coefficients a n and 
b n . This procedure is summarized in Table Al. 


# 
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Table A1 

Coefficients for the Stress Functions $ and y for the 

Isotropic Plate. 


Load 

a n 

b n 


rH 

II 

O 

O 

II 

c 

<T3 

b 2 = -1 


A -1 ■ 1 

a n = 0 

b 3 = -1 


A+l = 1 

* 

* 


a _2 s 1 

a n = 0 

b 4 - -1 


A +2 = 1 

• 

• 

ci 2 = 1 

• 

• 

b 4 = 3 

• 

• 


• 

• 

A - k * 1 

• 

• 

a n = 0 

• 

• 

V+2) * 

-i 

\k = 1 

\ * 1 

V*2) = 

k+1 


*See text 
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The equations which result from setting A + ^ = 1 and all other A n = 

0 are 


a 


V 


• b l° 


( A 1 5 ) 


0 = a^o 



(A 15 ) 


Thus, there are only two equations but three unknowns for A + j = 

1. A third equation is obtained when restrictions on the displacement 
field are imposed. The derivation of the displacement equation follows 
that of Chapter 3 with one modification. This modification is the 
result of the fact that the plate is assumed to be in a state of gen- 
eralized plane stress, while the pin was assumed to be in a state of 
plane strain. With this difference accounted for, the displacement 
equation becomes 


i v = 4> - - 4] 


(A17) 


where 



3-4v 
1-v * 


Recall that 



( A 1 8 ) 


(A19) 
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These Integrations yield a logarithmic term similar to that of eqs. 115 
and 116 of Chapter 4. The logarithmic terms are multivalued and must be 
accounted for in this analysis. Examining the displacements 
at 8 = 0 and 0 = 2n yields the following 

(u ' iv) e=2n = 2111 <*Vi ‘ V + (u ' iv) e=0 * (A20) 

If the displacements are to be single-valued, then the following must be 
true 


< a + b = 0 . (A21 ) 

Using equations A15, A16 and A21, the value for a^, b^, and b^ can be 
found. These are 


1 

a l = 

1 1+Oj 

( A22 ) 

K . 

1 

1 ‘ " 1+Kj ’ 

(A23) 

h 2 
b 3 " \+< l * 

(A24) 


With these results, the elastic response of the plate due to 
i k 6 

the e loads is completely defined. The analyses of Chapters 2 and 3 
can now be combined with this appendix to produce an algorithm capable 
of predicting the stresses in pin-loaded isotropic plates. The results 
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of this procedure will now be discussed. 

A standard case has been chosen to present the results for the 
isotropic case. Variations of the pin modulus, clearance, and the 
coefficient of friction will be made on the standard case. The para- 
meters adopted for the standard case are 

plate - Aluminum (E = 10 7 psi , v = 0.25) 

pin - Aluminum (E = 10 7 psi, v = 0.25) 

friction coefficient - 0.3 

loading direction - positive (l.e., p = -0.3) 

clearance - \ = 0.01 

rigid body pin displacement - 6 = 0.05 

All of these properties are felt to accurately represent actual values 
which may be encountered in typical aerospace applications. Variations 
on the standard case will now be presented. 

Figures A1 and A2 show the stresses at R = 1 for the standard case 
with variable pin modulii. Rigid, steel, and aluminum pins are com- 
pared. It is evident that pin modulus plays a small role In the results 
presented here. A rigid pin assumption appears to be well justified 
based on these findings. 

The stresses at R = 1 for the standard case with variable pin/plate 
clearance are presented in figs. A3 and A4. The changing contact angle 
is seen to alter the location of the peak hoop stress for these cases. 
This may affect the location of first damage in isotropic plates. 


0 (degrees) 


Fig, A1 Contact Stresses at R = 1 for the Standard Isotropic 
Case with Variable Pin Moduli i. 




6 (degrees) 


Fig. A2 Hoop Stresses at R = 1 for the Standard Isotropic 
Case with Variable Pin Moduli i. 
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Finally, figs. A5 and A6 show the stresses for the standard case 
with variable friction coefficients. In each case, the plate is assumed 
to be under monotonical iy increasing loads (i.e., ^ is negative). 
Friction is seen to alter the normal stress distribution and it causes 
the hoop stress to change from tensile to compressive at 0 = 0°. 

Other variations on the standard case will not be presented here 
due to their similarity with the results presented in Chapter 6. In 
fact the conclusions drawn for the orthotropic case appear to be valid 
for the isotropic case as well. 
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